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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

Space  Shuttle  launches  deposit  massive  quantities  of  hydrogen  chloride 
(HCi)  into  the  environment  (Reference  l).  At  present,  the  short-  and  long¬ 
term  environmental  impacts  of  the  HCI  releases  are  not  well  characterized, 
and  quantitative  monitoring  of  the  transport  of  the  HCI  Into  the  areas 
surrounding  the  launch  site  is  desirable.  Previous  models  and  measurements 
demonstrated  that  the  amount  of  HCI  deposited  on  the  ground  depends  on  the 
local  meteorological  conditions,  temperature,  humidity,  windspeed  and 
direction  (References  1-11).  Furthermore,  the  HCI  can  be  found  in  both  the 
vapor  and  aerosol  (hydrochloric  acid  droplets)  forms.  The  relative 
partitioning  between  the  vapor  and  aerosol  forms  depends  on  a  variety  of 
factors,  such  as,  the  HCI  and  water  vapor  partial  pressures,  temperature, 
and  the  concentration  and  composition  of  the  aerosol  nucleation  sites. 
Since  the  concentrations  of  the  HCI  in  the  vapor  and  aerosol  forms  may  be 
quite  different,  it  is  necessary  to  monitor  both. 

B.  BACKGROUND 

Several  HCI  monitors  have  been  previously  developed  (References 
12-16):  however,  none  are  capable  of  monitoring  both  Ihc  vape.'  and  aerosol 
HCI  content  in  real-time.  Spectral  Sciences,  Inc.  (SSI)  has  developed  a 
combined  HCI  vapor-aerosol  monitor  under  Phase  II  SBIR  (Small  Business 
Innovative  Research  program)  funding  from  the  Air  Force  Fngineering 
Services  Center  (AFESC),  Tyndall  AFB,  Florida.  The  basic  operation  of  the 
SSI/AFKSC  HCI  monitor  is  depicted  in  Figure  1,  and  a  summary  of  its  key 
operational  specifications  is  given  in  Table  I.  The  instrument  consists  of 
four  major  components;  (1)  an  air-sampling,  optical  measurement  unit,  (2) 
a  power  supply  and  analog  signal  output  unit.  (B)  an  analog- to-digital 
(a/d)  converter  and  timing  pulse  unit,  and  (4)  a  personal  computer  (PC). 
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Ficure  1.  Schematic  ot  the  Operation  of  the  SSl/AFESC  HCl  Monitor. 


TABLE  1.  SPECIFICATIONS 

FOR  THE  SSI/AFESC  HCl  MONITOR. 

.yiiumiiin  Detectable  HCl  Cone. 

100  ppb 

Dynamic  Rariije 

0.1  -  10^  ppm 

Absorption  Path  Lc;ni;th 

7  m  {variable  from  4-10  m) 

Detector 

PbSe  (cooled  to  -20  C) 

Sample  Averaging  Time 

30  s  (variable  from  O.l-lO^s) 

Sample  Cell  Refilling  Time 

30  s 

;  Electrical  Reejuirements * ^ ^ 

250  W/120  V  AC 

(a)  This  does  not  include  the  PC. 
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An  air  s.'..,ple  is  drawn  fhrouj^h  the  aerosol  evnaporai  or  winch  vapo'  :zis  tia' 
aerosols  and  releases  the  dissolved  HCl  Into  the  [jas  nhasi; .  By  turiinij;  tin- 
aerosol  heater  on  and  off  one  can  determine  both  the  total  HC 1  (aerosol  .ini 
vapor)  and  the  vapor  phase  HCl .  The  HCl  concentration  ineasiireneni  is 
performed  with  a  special  infrared  liiiht  source.  Data  acquisition  ari’! 
analysis  are  controlled  by  a  remotely  locaieci  PC.  A  1000  foot  siiqia!  caiiie 
allows  the  PC  to  be  Ificated  in  a  bunker  or  control  room  and  enatiies  ri-a  i 
time;  monitorinii  of  the  concentration  meiisiirenient . 


C .  APPROACH 


The  uevice  is  based  on  measuring  the  infrared  absorption  ol  ;jas -piiasc 
HCl  in  the  3.4  (itti  waveiericth  rejjion.  The  measuremc-nt  approacii  elite's  t  i  or 
othe-  infrared  devices  (References  14-lfi)  in  that  it  iiti  ii/i-s  a  nove  ■ 
ml  raced  iiqrit  source.  in  essence,  the  source  consists  ol  HCi  in  a  heated 
ceil  whicn  e.tiits  the  characteristic  HC!  line  spectrum  Tin;  pr’ii'a'v 
auvantaoe  of  usinp  a  i  ii^tit  source  spectrally  matclieu  to  the  aosorption 
spectrum  is  one  o!  species  selectivity.  T.his  enabl'S  quani  i  tat  :  v.' 
measurement  of  very  small  species  concentrations  in  liie  nresenre  ot 
strontily  absorbinii  species  whose  absorption  line  positions  are  spint  r;i  i  1  y 
iincorre  1  at  ea  with  the  source  line  emission. 

The  basic  measurement  tochnitjut;  is  summai-i zed  in  Fijjurc;  2.  TJie  .HCl 
lamp  emission  consists  ot  hi  gb- i  ntens  i  t  y ,  narrow  HCl  emission  nni's 
supi.T  imposed  on  a  1  ow- 1  ntens  i  ty  continuum.  The  continuum  emission  is  nue 
to  the  quartz  cell  cnnfininc  the  hot  HCl  i;as .  Typically,  for  the  HCl  ianips 
produced  to  date,  the  total  (integrated  over  wavelength)  HCl  line  eniission 
IS  appro.x  I  .ma  i  e  1  y  etpial  t<i  the  1r»tal  continuum  emission.  For  a  conventional 
infrared  source,  sucii  iis  a  continuous  black  body  c-niitter  (lieateo  wire),  t  iie 
total  emitted  energy  within  the  spectral  width  ol  the  HCl  lines  is  only 
around  1  percent  ol  the  total  black  body  continuum.  This  means  tliat  t  !ie 
use  ol  the  HCl  lamp  retluces  tlie  itite^-fcience  of  uni'orre  i  a  t  (.-d  atisorbers  iiv 
approx  1  mat  i;  ly  two  orde’-s  of  magnitude. 

The  HCl  .monitor  uses  a  single  pathlength,  gas  correlation,  spectra: 

(  ;  !  t  ei- ng  scheme.  The  iam{)  emission  is  passed  i  hrough  a  filter  wnee  i 
containing  small  quartz  ('eils  whie.'i  either  contain  an  HC  I  '  n  i  t  rngen  m’xture 
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( cit;s i incited  HCl  iibsorbinc)  or  just  nitrogen  (designated  nonubsorbing )  .  The 
HCl  concentration  in  the  absorbing  filter  is  adjusted  to  absorb  out  the  fICJ 
line  emission.  The  difference  between  the  total  integrated  signals  for  the 
two  filters  is  equal  to  the  total  HCl  source  line  emission.  This  quantity 
IS  (it.-termi  ned  when  no  HCl  is  present  in  the  sample  absorption  path  (White 
C(.‘il)  and  is  stored  in  the  PC  as  the  "zero"  calibration  constant.  When  HCl 
IS  introduced  into  the  sample  absorption  path  the  difference  signal  for  the 
two  filters  is  reduced.  The  difference  signal  can  be  quantitatively 
reJated  to  the  HCl  concentra;.ion  in  the  sample  path. 

The  use  of  a  correlated  light  source  in  conjunction  with  gas 
correlation  filtering  is  a  significant  improvement  over  the  use  of  a 
continuum  black-body  source  with  gas  correlation  filtering.  It  further 
reduces,  by  about  two  orders  of  magnitude,  the  sensitivity  of  the 
measurement  to  uncorrelated  absorbers.  For  example,  with  conventional  gas 
correlation  detection  it  would  be  difficult  to  detect  1  ppm  HCl  in  the 
presence  of  100  ppm  CH^  which  also  absorbs  strongly  in  the  same  spectral 
region  as  HCl.  However,  by  using  the  correlated  light  source  it  would  be 


DIFFERENCE 
S I GNAL 

(CAL  I  BRAT  I  ON  ZERO) 


DIFFERENCE 
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Figure  2.  Schematic  of  Basic  Measurement  Approach. 
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possible  to  detect  1  ppm  HC I  in  the  presence  of  up  to  10,000  p[)m  . 
Another  aovantajie  relates  to  the  required  accuracy  and  lony-term  stability 
of  the  zero  absorption  calibration  point.  For  conventional  gas  correlation 
nctection  the  zero  calibration  must  be  both  accurate  and  stable  to  better 
than  0.03  percent  of  the  absolute  signal  level  to  measure  an  HCl 
concentration  of  1  ppm.  However,  by  using  the  correlated  light  source  the 
accuracy  and  stability  of  the  calibration  need  be  only  good  to  better  than 
1  pert;ent  in  order  to  detect  3  ppm  of  HCl. 

The  next  section  presents  an  overview  of  the  instrument  subsystems  and 
the  calibration  approach.  More  detailed  discussions  on  the  actual 
operation  and  servicing  of  the  HCl  monitor  are  given  in  the  appendices: 
Appendix  A  -  Instrument  Performance  Model;  Appendix  B  -  HCl  Monitor 
Software;  Appendix  C  -  Optical  Train:  and  Appendix  I)  -  Electronic  Design 
for  the  HCl  Monitor. 
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SECTION  II 


INSTRUMENT  DESCRIPTION 


A.  OPTICAL  LAYOUT 

The  optical  layout  of  the  HCl  monitor  is  shown  in  Figure  3.  Light 
from  ttie  HCl  lamp  is  collected  by  a  calcium  fluoride  lens  (25  mm  diameter. 
50  mm  focal  length)  which  focuses  the  light  onto  the  filter  wheel.  A 
second  calcium  fluoride  lens  (25  ram  diameter,  38  mm  focal  length)  refocuses 
the  light  at  the  entrance  aperture  to  the  White  Cell  (Reference  17).  The 
White  Cell  mirrors  have  a  focal  length  of  25  cm,  and  the  effective  light 
collecting  diameter  of  the  rear  mirrors  is  50  cm.  The  White  Cell  mirrors 
were  machined  out  of  aluminum  blanks  which  were  plated  with  nickel  before 
polishing  and  overcoated  with  gold.  The  number  of  passes  in  the  White  Cell 
is  ati justed  to  28  for  a  total  pathlength  of  7  meters.  A  third  calcium 
fluoride  lens  (25  mm  diameter,  38  mm  focal  length)  refocuses  the  output  of 
the  White  Cell  onto  a  thermoelectrically  cooled  (two-stage)  PbSe  detector 
(3mm  X  3mm  square).  A  band-pass  filter  is  placed  in  front  of  the  detector 
to  isolate  the  most  strongly  absorbing  lines  of  HCl.  The  temperature  of 


Figure  3.  Schematic  of  HCl  Monitor  Optical  Layout. 


the  detector  is  controlled  actively,  using  a  feedback  circuit.  Tins 
circuit  uses  a  thermister  built  into  the  detector  package  to  sense  the 
detector  temperature  and  sends  a  correction  current  to  the  thermoelectric 
cooler,  which  is  proportional  to  the  difference  oi  the  actual  and  preset 
detector  temperatures. 

An  airtight  enclosure  is  placed  over  the  White  Cell  (see  Figure  4) 
which  allows  the  air  to  be  drawn  through  the  aerosol  evaporator  and  into 
the  sample  absorption  path.  The  sampling  pump  has  a  maximum  tiiroughput  of 
94  cc/s  and  a  power  requirement  of  15  watts.  Given  the  enclosure  volume  of 
2400  cc,  it  takes  about  26  seconds  to  replenish  the  sample  absorption 
volume.  The  replenishing  time  constant  is  the  fundamental  limit  to  the 
time  response  of  the  HCl  monitor  to  fast  changes  in  the  HCl  concentration. 
If  a  faster  response  is  desired,  a  higher  throughput  pump  can  be  used.  If 
one  is  only  Interested  in  the  vapor  phase  HCl,  then  the  White  Cell 
enclosure  can  be  removed  and  the  limiting  time  response  is  the  time  it 
takes  the  chopper  wheel  to  make  a  single  revolution  (about  0.05  seconds). 


Figure  4.  Picture  of  Optical  Components  With  (Right)  and  Without 
(Left)  Sample  Absorption  Path  Cover  and  Aerosol 
Evaporator . 
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HCl  LAMP 


3  . 


A  schematic  ot  the  HCL  liyiit  source  is  siifiwn  in  Figure  5.  A  (juariz 
tube  (5.0  cm  Lotui  by  1.8  cm  ammeter)  is  filled  with  an  HC]/nitroccn 
mixture.  Ttie  front  window  of  the  emission  cell  is  an  extremely  thin  quartz 
buotiie  (approximately  0.1  mm  thick)  developed  by  G.  Finkenbeiner  (Refer'cnce 
18).  Because  of  its  thinness  this  window  is  very  transmissive  (see  Fif'ur-f; 
G)  and  rias  a  very  small  emissivity.  The  thickness  ol'  the  rear  window,  1.5 
mm  tliick  quartz  flat,  is  not  critical  because  this  part  of  tlie  cell  is  not 
heated.  A  single  layer  of  insulating  material  (ceramic  fiber  paper)  is 
used  to  separate  tfie  cell  from  the  heating  wire  (26  gauge  nichrome).  T)ie 
purposf;  of  this  is  to  more  uriiformly  distribute  the  heat  and  background 
infi'ored  emission  from  the  locally  hotter  and  brighter  nichrome  wire. 
About  12  primary  nichrome  wire  "fingers"  are  spaced  uniformly  around  tLie 
ceil.  These  are  covered  by  a  layer  of  insulation  and  an  additional  six 
secondary  nichrome  "fingers"  are  used  to  compensate  for  radiative  and 
convective  losses  at  the  front  of  the  cell.  This  is  important  because  it 
prevents  tfie  Ht:  I  line  emission  from  becoming  self-absorbed  by  the  cooler 
gas  near  the  front  window,  thus,  reducing  the  sensitivity  of  the  system  to 
s.mail  concentrations  of  HCi  in  the  sample  absorption  path.  The  cell  is 
tiien  wrapEied  with  additional  Insulation  to  approximately  7  cm  diameter. 
Trie  insulation  is  baked  to  drive  out  ail  organic  binders  which  increases 
its  insulating  propertu.-s. 

A  thermocouEile  is  placed  close  to  the  nichrome  windings  and  is  used  to 
mainruin  t.he  cell  at  constant  temperature  over  a  long  time.  The  nominal 
operating  temperature  of  the  lamp  is  approximately  1000  K,  although  the 
exact  temperature  is  not  too  important  as  long  as  it  remains  constant.  Tlie 
I  ami)  consumes  about  12  watts  at  steady  state. 

Tin;  lamp  emission  spectrum  was  measured  with  a  1/2  meter  monochromator 
and  is  compared  in  Figure  7  to  model  predictions  based  on  the  known 
spectroscopic  and  radiative  transport  properties  of  HGi  (References  17-22). 
From  an  analysis  of  the  lamp  e.mission  s[)ecrrum  it  was  determined  that  the 
elli.-ctivc  emissivity  of  the  (luartz  was  0.02  and  that  the  integrated 
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.ToJerular  lini;  anti  coniiniium  emissions  over  the  measurement  band  pass  were 
nearly  equal.  The  measurement  band  pass,  chosen  to  select  the  lamp 
emission  lines  which  are  most  strongly  absorbed  by  the  much  cold(;r 
atmosphet ic  HCl  .  is  shown  in  Finure  8. 


ALUMINUM 

HOUSING 

MACHINED 
CERAMIC  SUPPORT 

SEALED  OUARTZ 
CELL  CONTAINING 
HCl/N^  MIXTURE 


THERMOCOUPLE 


Fiyure  5.  Schematic  of  the  HCl  Emission  Lamp. 


EMISSION  CELL  WINDOWS 


WAVELENGTH  (Vn) 

Fijiure  6.  Transmission  of  Quartz  Windows. 
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Figure  7.  Comparison  of  Observed  and  Predicted  HCl  Lamp  Spectrum. 


3.2  3.M 


WAVELENGTH  (MICRONS) 


Figure  8.  Transmission  Spectrum  of  the  Band-Pass  Filter. 


c. 


FILTER  WHEEL 


The  filter  wheel  and  its  drive  system  are  shown  in  Ficure  ‘J .  It 
contains  four  quartz  cel's  (25  mm  diameter,  8  mm  thick)  of  which  two  are 
filled  with  ary  air  and  two  are  filled  with  an  HCl/nitrojjen  mixture.  Ttie 
conc:entrat  i  on  of  HCl  in  the  filter  cell  is  adjusted  so  that  most  of  tlie  HCl 
line  emission  from  the  lamp  is  attenuated.  Predictions  of  the  attenuation 
of  the  lamp  emission  for  several  HCl  concentrations  in  the  filter  cell  are 
displayed  in  Figure  10.  The  current  version  of  the  prototype  HCl  monitor 
uses  an  HCl  mole  fraction  of  0.11.  The  two  kinds  of  cells  are  arrarigeu  so 
that  like  cells  are  180  degrees  opposed.  It  is  not  necessary  to  use  a 
reference  pulse  to  distinguish  among  the  cells  or  to  indicate  the  overall 
phase  of  revolution  of  the  wheel.  This  is  accomplished  by  the  PC  signal 
processing  software  which  recognizes  that  the  signal  alternates  between 
just  two  different  maxima  with  the  largest  maximum  always  corresponding  to 
the  pure  nitrogen  filter  cell  (largest  transmission).  The  filter  wheel 
rotates  at  20  revolutions  a  second;  this  rate  is  held  constant  by  a  control 
circuit.  This  circuit  uses  a  detection  output  to  sense  the  rotation  rate 
and  provides  a  correction  voltage  to  the  motor  which  is  proportional  to  the 
aifference  of  the  actual  and  preset  rotation  rates. 


Figure  9.  Schematic  of  Filter  Wheel. 
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CD 


WnVELENGTH  (Uf^) 


Figure  10.  Predicted  Attenuation  of  HCl  Lamp  Emission  for  Several 
HCi  Concentrations  in  the  Filter  Cell.  Filter  Cell 
Conditions:  Temperature  (295  K),  Total  Pressure  (0.9 

atm),  and  Pathlength  (1  cm). 

D.  SAMPLE  CELL  TRANSMISSION 


The  attenuation  of  the  HCl  lamp  for  various  concentrations  of  HCl  in 
the  sample  ceil  was  experimentally  determined  using  the  1/2  m  monochromator 
at  high  spectral  resolution  so  that  the  lines  of  the  individual  chlorine 
isotopes  (H'^^Cl  and  H^^Cl)  were  resolved.  This  is  shown  in  Figure  11  for 
one  of  the  absorption  lines  within  the  measurement  band  pass.  A  1  cm  path 
length  absorption  cell  was  used  for  these  spectral  scans.  The  equivalent 
concentration  for  the  7  meter  pathlength  of  the  actual  instrument  is  equal 
to  1/700' t)i  of  tlie  concentration  in  the  1  cm  cell.  Model  predictions,  also 
shown  in  Figure  11,  are  in  good  agreement  with  the  experimental 
measurements  and  confirm  that  even  small  concentrations  of  HCl  in  the 
sample  path  can  appreciably  attenuate  the  lamp  output.  About  a  1  percent 
reuuction  in  source  intensity  occurs  for  every  ppm  of  HCl  in  the  absorption 

f 

path.  Ttiis  can  be  compared  to  Instruments  based  on  a  continuum  infrared 
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source,  such  as  the  Lawrence  Livermore  laboratory  Instrument  (Reference 
16),  where  there  Is  approximately  a  0.04  percent  reduction  in  source 
intensity  for  every  ppm  of  HCl  in  the  absorption  path  (this  assumes  a  7- 
meter  total  path) . 


REFERENCE  CELL  RBSORPTION  CELL 


3.370  3.375  3.380  3.370  3.375  3.380 


WRVELENGTH  (MICRONS)  URVELENGTH  (MICRONS) 


Ficure  11.  Comparison  of  Observed  and  Predicted  Attenuation  of  A 
Single  Pair  of  Isotopic  HCl  Emission  Lines.  Lamp 

Conditions;  Temperature  (1100  K),  Total  Pressure  1.01 
atm),  Pathlength  (4.5  cm),  HCl  Mole  Fraction  (0.1). 
Quartz  Emisslvity  (0.02).  Equivalent  White  Cell 
Conditions : Temperature  (295  K),  Pressure  (0.92  atm).  Path 
Length  (7  ra),  and  HCl  Concentration  (32.9  ppm). 
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AEROSOL  EVAPORATOR 


£ . 


Tru;  aerosol  evaporator  is  displayed  in  Figure  12.  Air  is  drawn 
through  the  resistiveJy  heated  porous  carbon  plug  (2.54  cm  o.d,  1.91  cm 
i . d .  )  .  The  number  of  times  an  aerosol  particle  encounters  the  internal 
siirtaces  of  ttie  carbon  plug  as  it  passes  through  the  matrix  is  aetermined 
by  the  tiiickness  and  porosity  of  the  carbon.  The  total  power  dissipated  by 
th(f  lu;ater  is  about  22  watts  (4.7  volts  through  1.0  ohm). 

The  steady-state  temperature  of  the  gas  exiting  the  heater  is  IIS^^C. 
Ttiis  t(;mperature  is  attained  in  less  than  10  seconds  after  switching  on  the 
room  temperature  evaporator.  However,  there  is  a  long  decay  time  down  to 
room  temperature  after  the  evaporator  is  switched  off.  It  decays  to  54°C 
in  20  seconds  and  takes  an  additional  6  minutes  to  decay  to  within  5°C  of 
!’oom  temperature.  This  long  decay  time  is  due  to  the  solid  carbon  base, 
the  metal  end  caps  and  the  stainless  tubing  which  get  quite  hot  when  the 
evaporator  is  on  and  cool  slowly  by  transferring  heat  to  the  incoming  air 
afrer  the  evaporator  is  switched  off. 
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F. 


DATA  ACQUISITION  AND  SIGNAL  PROCESSING 


The  PbSe  detector  is  operated  in  a  standard  ac  coupled  photo- 
conductive  mode.  The  detector  signal  is  amplified  and  buffered  to  provide 
a  low-impedance  sicnal  output,  which  is  proportional  to  the  incident  liifht 
intensity.  To  facilitate  transmission  of  the  ac  signal  through  the 
1000-foot  shielded  coaxial  cable,  a  6-voit  dc  pedestal  is  added  to  the 
detector  output.  The  transmitted  signal  is  again  ac  coupled  to  remove  the 
6-volt  bias  and  the  resulting  analog  signal  is  shown  In  Figure  13.  A 
sequence  of  timing  pulses  is  generated  to  trigger  the  a/d  converter.  The 
timing  pulses  arc  generated  using  a  circuit  which  outputs  a  pair  of  voltage 
spikes  by  differentiating  the  detector  signal.  Both  the  time  separation  of 
the  two  spikes  and  as  the  time  delay  of  the  first  spike,  relative  to  the 
signal  leading  edge,  can  be  Independently  controlled.  This  enables  the  two 
spikes  to  be  precisely  positioned  on  the  signal  curve,  one  at  the  maximum, 
and  the  other  at  the  minimum.  Changes  in  the  rotation  rate  of  the  filter 
wheel  will  cause  the  timing  pulses  to  occur  at  positions  offset  from  the 
maximum  and  minimum.  This  is  why  the  rotation  rate  is  controlled. 
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Figure  13.  Schematic  of  Analog  Output  and  A/..  Triggering  Signals. 
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Tin;  siijiujjs  ;iri;  ac(iuin;d.  proc(;ssod  in  rnai  -  tine,  .'irui 

curiv(.‘rrea  to  an  HCl  concentration  using  FOKTRAN- based  software  running  on 

an  [BM -coiigia  t  i  b  lo  PC.  A  difference  is  taken  between  each  digitized  iraxi'nurT! 
and  tile  average  of  its  adjacent  digitized  baseline  points.  This  produces 
two  signals,  one  designateu  which  correspontis  to  tlie  nitrogen  filter  ane 
IS  the  most  sensitive  to  HCl  absorption  in  the  sam^ile  cell,  and  the  ottier 

designateu  wliicii  corresponds  to  the  HCl  filter  and  is  essentially 

nonabsorbeu  by  the  HCl  in  the  absorption  path.  A  ratio  of  these  signals  is 
co.niputed  , 


w.'iere  Liie  difference  signai  is  noriraiized  to  the  nonabsorbing  signal.  This 
normalization  mitigai.es  tiie  effects  of  signal  drifts  caused  by  f  1  uc  t  ua  t  i  oris 
'.  n  tin;  taitp  oiitpul  ane  accii.mu  Jat  ion  o!  surface  cont  am  i  na  t  i  on  on  ti;e  White 
(c'li  mirrors.  An  oveTa!.'  syste.m  calibration  constant,  ,  is  determined  by 
.measuring  trie  above  rat'ci  witliout  HCl  in  tlio  sample  ciiamber.  The 
concentration  of  iiCi  is  then  related  to  a  normalized  ral  lo , 


wncre  this  quantity  is  equal  to  unity  for  zero  HCJ  concentration  and 
approaciies  zero  for  large  HCl  concentrations.  This  ratio  is  actually  the 
average  ot  the  ratios  I'or  a  number  of  wheel  rotations.  The  number  of 
roiat.ons  chosen  to  be  av(;raged  is  artiitrary  and  is  selected  by  the  user. 
Ty[)icaily,  400  wheel  rotations  are  averaged  wiiich  takes  about  30  seconds 
and  corresponds  to  the  replenishing  time  of  tiie  sample  absorption  volume. 

In  del  (.'rmining  the  HCi  concentration  from  the  .measured  intensity  ratio 
il  is  conven  i  en L  to  work  with  the  ipiantity 


wiiicii  IS  proportional  to  the  coneentrat  ion .  In  theory  (Reference  22) 
cnnceM I la t  1  on  can  be  related  to  this  quanlity  via 


A,X  + 
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where  C  is  concentration,  ami  ar-e  co<!nicionts  (n;t(!rni  i  iieci  t)y 

experimental  calibration.  Hichcr -order  terms  can  be  added  if  deemed 
necessary.  The  effect  of  temperature  on  the  concentration  determ i na i i on 
must  also  be  inchuieii.  Atmospheric  temperatures  can  typically  ranine  fr(i"i 
arouno  0^^  C  ourtng  the*  winter  to  as  high  as  50‘^’C  on  a  iiot,  sunny  nay  in  the- 
nesert  Also.  the  aerosol  evaporator  nischartjes  j^as  into  the  sanntie 
atisorption  reiiinn  at  around  The  effect  of  temperature  can  ue 

empirically  included  in  the  expansion  coefficients  ttiroU);n 

(b) 

evainated  at  the  relerence  tempera  i  tire ,  1',^ . 

detei'mined  from  model  predictions  of  the 

G.  CALIBRATION' 

The  eaiihration  coefficients  wtire  determined  at  room  temperature, 
us  I  nj;  Known  mixtures  of  HC 1  in  nitro);en  wnich  were  pre[iared  in  a  1  cm 

p<it  ii  leiiijl  1:  ceil.  This  1  cm  cell  was  then  positioned  at  the  r-x  i  t  of  the 
sample  aosor'pt  id.n  cell  arul  the  intensity  ratio  (as  in  Eifnation  (2))  was 
dcterm  t  ni‘d .  The  equivalent  concentration  in  th.e  sam[)le  path  whieli  would 
result  in  this  same  intensity  ratio  would  be  !/70(i'th  (the  ratio  of  path 

lencths)  or  t  ne  conci'nt  rat  ion  of  the  !  I'm  ceil.  By  varyinj;  the 

concf-nt  ra  t  1  on  of  1  in  the  1  cm  ct;  t  I  a  larjo-  ranjte  of  eiiuivaleiit  sanpie 
cell  concentrations  were  sim.ulaied.  The  HC  1  coneent  rat  ion  tn  t  lie  1  cm  ceil 
was  uetermined  by  measuring;  the  speet  rally  resolved  1  r-ansm  i  t  i  ance  fur  a 
niack  oouy  lu;l"  source,  t  lien  tisinn  a  synthetic  sriectrai  model  to  match  trie 
measured  traiisrr  i  t  vutuie .  An  exar.ipie  of  tin;  result  ut  this  procedure  is 
s.'idwn  Ml  ■■'.jjure  ;  ■)  .  Tin.-  cMiibration  curve  was  defi>rm:niMi  in  two  stiqis. 
First  tlie  nriuts  to  the  theoretical  .modc-l  were  varied  from  tlieir  noiinai 

vaiues  to  ijive  a  yood  tit  to  the  calibration  points.  Tiuin  the  t  .heore  t  ica  1 
curve,  wnich  was  determined  at  a  much  finer  firid  spacini^  than  the 

eaiihration  points,  was  tit  to  a  sixth-order  polynomial  expansion  mi  X,  as 
in  ruination  (2).  The  result  inij  calibration  cur-ve  is  shown  :n  Fiijure  15, 


'  I  o 

■  o 

w.’ii'rt-  C,|  ’s  I  ne  concent  I'a  I  3  on 
and  a  is  a  s'-ai'icj  exponent 
t  e':qie>-ai  ui  e  var  i  at  i  on  . 
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The  much  larger  error  bars  on  the  two  lowest  concentration  calibration 
points  reflect  the  fact  that  the  concentration  for  these  points  could  not 
be  directly  measured.  Because  their  transmi ttances  were  too  close  to  unity 
to  yield  a  measurable  absorption  loss,  the  concentrations  were  estimated  by 
using  the  relative  dilution  factor  established  for  the  other  higher 
cmicentra t ! on  points.  The  overall  accuracy  of  the  calibration  is  estimated 
to  be  approximately  ±  20  percent. 

The  temperature  correction  exponent  was  determined  by  using  a  first- 
principles  model  of  the  instrument  detector  output.  This  model  was 
verified  to  give  good  agreement  with  the  experimentally  determined  room 
temperature  calibration  curve  and,  therefore,  can  be  used  with  confidence 
to  extrapolate  the  calibration  curve  to  higher  temperatures.  Examples  of 
the  temperature-dependent  calibration  curves  over  the  full  temperature 
range  of  interest  are  shown  in  Figure  16.  The  temperature-dependence  is 
relatively  small  and  is  mostly  due  to  the  decrease  in  HCl  density  at 
constant  pressure  as  the  temperature  is  increased. 
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Figure  14.  Comparison  of  Observed  and  Predicted  Transmission  of  a 
Hlackbody  Source  Through  an  Equivalent  Sample  I'ath  HCl 
Concentration  of  147  ppm. 
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HCI  CONCENTRATION  (PPM)  HC I  CONCENTRATION  (PPM) 


An  initial  test  of  the  instrument  was  performed  usinu  precalibrated 
mixtures  of  HCl  in  nitrotjen  which  were  supplieu  by  Matheson.  The  results 
are  shown  in  Ficure  17.  The  nominal  concentrations  for  the  two  mixtures 
were  rl  and  10  ppm.  We  have  found  with  previous  HCl  uas  mixtures  that  ttie 
concentration  can  change  substantially  and  unpredictably  over  lonu  perious 
of  time.  The  particular  mixtures  used  were  about  8  months  old  so  that 
exact  aj'reenient  with  the  nominal  concentrations  was  not  expected  and  the 
actual  re.sults  are  quite  reasonable.  The  initial  measurement  points  up  to 
approximately  3  minutes  were  room  air  samples  and  indicate  the  measurement 
noise  for  a  15-second  measurement  averaue.  When  the  61  ppm  HCl  sample  was 
initially  turned  on  there  was  a  larye  concentration  spike.  This  was  due  to 
outqassing  from  the  regulator  which  had  previously  been  used  on  a  much 
hiijher  HCl  concentration  mixture.  The  concentration  decayed  rapidly  when 
the  HCl  sample  was  turned  off  at  about  22  minutes.  The  effective 
evacuation  time  for  the  sample  absorption  cell  is  around  1  minute  which  is 
consistent  with  the  pump  throughput.  The  10  ppm  HCl  mixture  was  turned  on 
at  27  minutes  and  an  initial  concentration  spike  was  again  observed  due  to 
regulator  outgassing. 


Matheson  HCI/N2  Mixtures 
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Figure  17.  .Measurement  of  Precalibrated  Matheson  HCI/N2  .Mixtures. 
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H.  OPERATION  AND  MAINTENANCE 


Operation  and  maintenance  instructions  for  tlie  HCI  monitor  arc;  brii-fly 
ucscribeo  in  this  section.  Oetailed  instructions  are  civen  in  Appendic(;s 
A-D.  Tilt;  loLJowinjj  step-by-step  directions  are  for  S(!t  up  and  operation: 

•  Ttie  wooden  stiippinc  crate  containinK  tin;  air 

samplintj/optical  m(;asurement  and  power  supply  modules  js 
placed  in  the  area  to  be  sampled. 

•  After  removing  tlie  crate  top  the  aerosol  heater  and  heavy 
duty  extension  cord  are  removed. 

•  The  extension  tube  for  tin;  aerosol  heater  is  attached  to 
the  air  intake  port  on  top  of  the  optical  moduli;. 

•  The  extension  cord  is  connected  to  the  crate  fan  and 
power  supply  power  cords. 

•  Turn  the  switch  for  the  main  power,  located  on  the  power 

supply  front  panel,  to  the  ON  position.  A  ijreen 

indicator  li}jht  will  turn  on. 

•  Turn  the  switch  for  the  aerosol  heater,  also  located  on 
the  front  panel,  to  ON.  A  red  indicator  light  will  turn 
on.  The  aerosol  heater  and  pump  can  be  controilod 
remotejy  from  the  computer  terminal. 

•  The  coaxial  signal  cable  is  attached  to  the  front  panel 

of  the  power  supply  module. 

•  The  crate  top  may  now  be  replaced.  Pull  the  two  wire 

cord  for  the  aerosol  heater  through  the  opening  in  tin; 
crate  top. 

•  The  aerosol  heater  is  attached  to  the  extension  tube. 

This  is  a  Swage  lok  fitting  and  should  not  be  excessively 
t ightened , 

•  The  (;]ectrical  connection  to  the  aerosol  heater  is  made 

with  the  polarized  two  pin  connector  pulled  through  the 
hole  in  the  crate  top. 


This  completes  tlie  set  up  of  the  sampling  unit.  The  control  unit  may 
be  remotely  located  from  this  module.  its  set  up  is  as  follows: 

•  Tlie  other  end  of  the  coaxial  cable  is  connected  to  the 
remotely  located  Ph . 
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rtii;  !iCl  priMirnni  (iisk(;tte  is  loadert  into  either  disk  drive 
ana  tfie  cuiuputcr  is  then  turned  on. 


•  After  the  eomputer  operating  system  is  loaded,  it  takes 
about  a  minute,  a  prompt,:,  will  appear  on  the  screen. 

•  Ttie  HCl  monitor  is  then  activated  by  typing  "hcl"(  ignore 
ciuotes)  after  the  prompt. 

•  A  series  of  questions  concerning  date  file  names,  data 
avt;ruging  and  measurement  time  period  are  then  answered 
by  the  user. 

•  Tile  resuits  of  tin;  measurement  are  displayed  in  real  time 
on  tiie  screen  in  the  form  of  a  graph  of  concentration,  in 
units  of  ppm,  versus  time. 

•  The  measured  data  is  also  placed  in  a  user  designated 
file  which  can  he  written  to  a  diskette  for  permanent 
storage . 


Tile  only  routine  maintenance  required  is  to  check  and,  if  necessary, 
to  cfean  the  White  cell  mirrors  after  each  field  measurement  program. 
Specific  instructions  for  doing  this  arc  given  in  Appendix  C.  Procedures 
for  changing  tiie  calibration  constants,  replacing  the  White  cell  mirrors, 
realigning  the  optical  system,  and  adjusting  the  electronic  control 
circuits  are  fully  detailed  in  Appendices  A-D. 
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SECTION  HI 


CONCLUSIONS 

SSI  has  successfully  developed  and  tested  a  novel  HCl  monitor  capable 
of  detecting  HCl  in  both  vapor  and  aerosol  phases.  An  innovative  HCl 
molecular  line  lamp  is  used  as  the  light  source  enabling  selective  and 
sensitive  uetection  of  HCl.  A  lower  detection  limit  of  0.1  ppm  HCL  has 
been  demonstrated  with  the  prototype  HCl  monitor.  The  monitor  has  been 
uelivered  to  AFESC  for  use  in  an  ongoing  program  of  monitoring  solid 
propellant  firings  at  several  Air  Force  test  and  launch  facilities. 
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APPENDIX  A 


INSTRUMENT  PERFORMANCE  MODEL 


The  detected  signal  can  be  expressed  by 


I  =  rcioTf  (o)Tg{u)T(,(cj){l-T^(o)T„(t>)  }N{cj.T2) 


(A.l) 


where  <j(cm~^)  is  wavenumber,  t  is  transmission  for  the  band-pass  filter, 
subscript  "f,"  for  the  sample  cell,  subscript  "s,"  for  the  gas  filter  cell, 
subscript  "c,"  for  the  emission  lamp,  subscript  ”1,"  and  for  the  quartz 
cell  bubble  window,  subscript  "w,"  and  N(o,T2 ) (W/sr/cm^/cm"^ )  is  the 
black-bocty  function  evaluated  at  the  lamp  temperature,  .  There  are 
frequency  independent  multiplicative  factors  such  as  the  collection  solid 
angle  which  are  required  to  convert  the  quantity  in  Equation  (A.l)  into  a 
quantitative  signal  level.  However,  these  factors  cancel  out  when  signal 
ratios  are  used  to  determine  the  concentration  of  the  absorbing  species  in 
the  sample  cell.  The  black-body  function  Is  computed  from 


N(u.Ti) 


'ir{exp(c2CJ/T2  )-l ) 


(A. 2) 


where  Cj ( 3 . 7405x10'^^ )  and  Cgd -43879)  are  the  first  and  second  radiation 
constants . 

In  general,  the  transmission  function  for  a  gaseous  absorber  can  be 
determined  from 


T  =  exp(-cplk(u,T) )  ,  (A. 3) 

where  c  is  the  mole  fraction  of  the  absorbing  gas,  p  (atm)  is  the  total 
pressure,  1  (cm)  is  the  path  length,  and  k  (cm"^  atm“^)  is  the  frequency 
ard  temperature-dependent  spectral  absorption  coefficient.  The  absorption 
cojfficient  for  a  single  molecular  line  is  given  by 
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kj(u.T) 


(tJ-Oj)2+Tj2 


(A. 4) 


where  Sj(cm“2  atin“^ )  is  the  Integrated  line  strength  for  the  jth  line. 
Tj(cm~^)  is  the  line  width,  and  (jj(cm~^)  is  the  line  center  location.  The 
total  absorption  coefficient  is  found  by  summing  over  all  the  kj's  that 
contribute  to  the  filter  band  pass  region,  and  this  is  the  k(u)  that  is 
used  in  Equation  (A. 3). 

The  line  width  depends  on  both  temperature  and  pressure  and  is 
calculated  from 


Tj  =  rj^p^T^/T 


vO  O 


(A.5) 


where  is  the  foreign  gas  broadening  line  width,  is  the  self¬ 
broadening  line  width,  and  Tq(K)  is  a  reference  temperature  (taken  to  be  at 
273  K  for  the  HCl  lines).  For  HCl  the  line  broadening  coefficients  depend 
on  the  particular  rotational  level  associated  with  the  transition  and  can 
be  approximated  by 


=  0. 097-7. 55xl0"2j  ,  and  U.6) 

t9  =  0.15-0.015j  ,  (A. 7) 

J  ® 


where  j  is  the  rotational  quantum  number  of  the  lower  level  Involved  in  the 

transition,  is  set  equal  to  0.014  cm~^  for  values  of  J  greater  than  11, 

and  is  set  equal  to  zero  for  j  greater  than  10. 

J  S 

The  line  strength  depends  on  the  temperature  and  is  computed  from 


“mSo<'">fm 

(Joqj.(T)qv(T) 


expf-CgByCm-l  )m)/T)  {l-exp{-C2to^/T} 


(A. 8) 


where  is  the  total  integrated  band  strength,  m  =  j+1  in  the  R  branch  and 
m  =  -j  in  the  P  branch,  <0q  is  the  band  origin,  (m)  denotes  the  absolute 
value  of  m,  0^,  is  the  transition  frequency  of  the  m'th  line,  fjj,  is  the 
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Herman-Wailis  intensity  correction  factor,  q^.  is  the  rotational  partition 
function,  is  the  vibrational  partition  function,  and  Bq(10.6  cm”^  for 
HCl)  is  the  ground  state  rotational  constant.  The  partition  functions  are 
given  by 


qr(T)  =  — and 

C2«o 


l-exp(-C2CJo/T) 

The  Herman-Wallis  factor  is  given  by 


(A. 9) 

(A. 10) 


f^  =  1. -0.025m  .  (A. 11) 

The  transition  frequencies  are  determined  from 

=  Oq  +  a^m  -  32™^  -  a3m^  ,  (A. 12) 

where  the  expansion  coefficients(a2 ,82.33)  and  the  band  center  location, 
Oq,  are  different  for  each  isotope  of  HCl.  These  values  are  tabulated 
below. 


different 

for  each 

isotope 

Of  HCl. 

These  values 

“0 

^1 

^2 

^3 

H^^Ci 

2886.04 

20.5725 

0.3072 

2.13x10"^ 

2883.92 

20.5425 

0.3065 

2.13x10"^ 

The  integrated  line  strength  for  HCl  is  155  cm"^  atm“^  however  before  using 
it  in  Equation  (A. 8)  it  must  be  adjusted  to  account  for  the  isotopic 
abundances.  Thus  for  H^'^Cl  use  Sq  =  38.8  and  for  H^^Cl  use  Sq  =  116.2. 

The  transmission  of  the  quartz  bubble  window  is  given  by 


=  epx(-x„) 


(A. 13) 


where  x^^  is  tfie  optical  opacity  of  the  window(  typical  ly  x^=0.02). 
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The  Detected  signal  given  in  Equation  (A.l)  is  dKicijnined  tor  tour 

sets  of  conditions:  (1)  no  HCl  in  the  sample  absorption  path  and  the 

nitrogen  containing  filter  cell  in  the  beam  path,  designated  1°  (2)  no  HCl 

<1 

in  the  sample  path  and  the  HCl  filter  cell  in  the  beam  path,  1°^^,  (3)  HCl 
in  the  sample  path  and  the  nitrogen  filter  cell  in  the  beam  path,  and 

(4)  HCl  in  the  sample  path  and  the  HCl  filter  cell  in  the  beam  path, 

The  first  two  signals  are  used  to  compute  the  zero  calibration  constant  via 


=  }/l“ 

’a  na  na 


(A. 14) 


In  general,  this  calibration  constant  depends  on  the  temperatures,  HCl 
concentrations,  total  pressures,  and  path  lengths  of  both  the  HCl  lamp  and 
the  filter  cells.  However,  for  a  given  instrument  all  these  factors  are 
held  constant  except  for  the  temperature  of  the  filter  cells  whose 
temperature  is  determined  by  convective  cooling  of  the  optics  can. 
Variations  in  the  ambient  temperature  will  then  lead  to  variations  in  the 
filter  cell  temperature.  This  is  why  the  zero  calibration  constant  must  be 
determined  in  real  time,  just  prior  to  an  HCl  release.  A  10‘^C  temperature 
variation  will  lead  to  approximately  a  1  percent  shift  in  the  calibration 
constant  and  would  produce  about  a  0.7  ppm  HCl  concentration  offset  in  the 
calibration  curve. 

The  HCl  concentration  in  the  sample  absorption  path  is  related  to  a 
measured  sig  al  ratio,  modeled  by 

R  =  —  ,  where  (A. 15) 

^o 


r 


^a~^na 


na 


(A. 16) 


In  constructing  a  calibration  curve,  it  is  convenient  to  work  with  the 
quantity 


X  =  -InR 


(A. 17) 
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which  is  proportional  to  the  HCl  concentration.  A  theoretical  calibration 
curve  can  be  constructed  simply  by  computing  X  as  a  function  of  the  HCl 
concentration.  C,  in  the  sample  absorption  path.  The  resulting  curve  can 
be  fit  with  a  simple  polynomial  expansion 


C  =  AjX  +  A2X^  +  A3X^  +  ■ •  • 


(A. 18) 


where  the  number  of  expansion  terms  depends  on  the  desired  accuracy  of  the 
fit  (typically  up  to  X^  is  required  for  the  current  application).  The 
calibration  curve  is  temperature  dependent  and  this  can  be  simply 
incorporated  into  the  via 


C  = 


(A. 19) 


where  is  the  concentration  evaluated  at  the  i  erence  temperature  T^  and 
a  is  a  scaling  exponent  determined  from  model  predictions  of  the 
temperature  variation.  For  this  application  a  =  3/4  was  found  to  yield 
good  approximation  to  the  temperature  dependence. 
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APPENDIX  B 


HCL  MONITOR  SOFTWARE 

The  HC I  monitor  is  a  m  i  rrocomputer- interlacorl  instrument  wriose 
sottware  is  oosoribea  in  tins  section.  As  was  merit  loneu  prirviousiy  in  riie 
eiecrronics  section,  two  siiiriais  •.■nter  the  IJATA  TKA.VSLATfO.VS  anaiosi  to 
(iijiilal  board.  Tite  first  signal  is  the  airpiifieu  output  of  tire  PbSe 
(U.-tector  Figure  B.la.  The  sircond  signal  is  a  clock  pulse  derived  from  the 
first  signal  which  tells  th«;  A/I)  boaro  wnen  to  read  Figure  3.1b.  Two  files 
are  supplieo  with  the  disk,  'hci.t.'xe"  and  "calib".  The  later  file  is  the  J 
ppm  HCl  calibration  constant  read  by  hcl.exe.  This  constant  is  upuateo 
each  time  a  calibration  run  is  performed.  The  former  file  is  tlie  actual 
HFl  monitor  data  collection  calibration,  and  analysis  program.  A  flowchart 
for  this  program  is  shown  in  Figure  B.2. 

Upon  entering  the  program  a  aisplay  ivlentifying  the  program  is 
flashed  to  t.he  screen  for  3  seconds.  Then  the  user  is  asked  whether  or  not 
a  calibration  run  is  to  be  performed  and  for  how  long.  11  is  ativlsable  to 
perlorm  an  instrument  calibration  run  for  at  least  1/2  hour  (after  tin; 
instrunK.'nt  warmup  period)  before  collecting  data.  This  procedtirt;  ensures 
that  the  a[)propriate  U  ppm  HCl  caiibrtition  constant  will  be  iisktn  in 
(letermining  the  HCl  concentrations  to  be  later  rnttasured.  The  calibration 
procedure  should  be  performed  while  t h-»re  is  no  HCl  present.  If  there  is 
sorni!  small  (a  few  ppm)  HC  1  concentriit  i  on  present  during  ea  It  brat  ion  tht'.n 
the  subseipient  measuring  period  may  report  relative  KC 1  concentration  ar.u 
not  absolute  concentration.  The  ca i i tiraf ion  constant  is  stored  in  the  file 
(jaiib.  "  at  the  enii  of  the  procedure  fur  later  recall  during  a  measurement 
rtin . 

Foi  lowing  the  caltbriition  of  the  instrument,  a  measurement  ritn  can 
be  :ier  f  or:;a;d .  After  selecting  tfie  number  of  wheei  turns  to  averagt;,  t.he 
user  is  [irompted  fui-  the  name  of  the  storage  file  and  the  duration  of  the 
run.  it.  IS  suggesled  that  data  files  be  stored  on  a  separate  fioriny  to 
protect  against  exceeding  disk  sptice  on  the  program  uisk  and  luince  losing 
rue  (lata.  The  program  opens  the  data  file  and  stores  some  user  suppiieu 
.header  informatjon  tit  it.  Then  a  graphic  display  of  HCl  conct'nt  rat  i  on 


vc'siKs  i;niL*  is  pui  mi  I  lit;  surmeii.  The  pfoi^r-ain  iiala  feuuc  t  i  mt  anu  ariaiysi:, 
omars  !  ii  rfie  suDroutira.'  "anaiOI."  h'ssent  i  a !  i  y  this  rourine  takes  the 
n  1  !  t  ei'eace  aetweeri  the  lai't^e  peaks  anu  the  small  peaks  ot  fii^ure  B.l  a  thea 
('.:\’!i;es  iiy  the  small  [leaks.  These;  are  then  averaj^en  over  the  wheel  tarns 
.iini  reiMirteU  bacK  to  tha;  main  [irtiyram  as  ttn;  nejpitive  luijar  i  t.han  at  this 
ip.aai  ity.  '."in;  mairt  [irojjr'ai;:  tfien  reiains  this  Loi^ar ;  t  hm  to  the  hh  1 
enaeeat  ra  t  ;  on  usitiji  a  caitoratjon  curve.  Durinji  the  actju  i  s ;  t  ;  mi  arui 
anaiyss  the  (t  iitie,  coaeent  ra  i  i  mi )  paints  are  (iis[)iaye(i  on  the  scr-eea  ana 
oaiiiat  to  the  data  file.  At  the  end  of  acnuisition  the  file  is  cioseti  and 
1  tie  tirofjrant  sto[)s.  At  any  Linn;  duriny  the  data  actjuisition  tirocess  tfie 
pro;;raiii  ettn  bi;  stop[)ed  and  the  file  saved  by  pressing  the  "escape"  key. 

A  projiram  list  me  follows  in  Section  B.l.  The  code  is  implemented  in 
;i  mi.xture  of  .PORTRAN',  ,  and  ASSi'MBhKR.  The  FORTRAN  modules  were  (:om[)il('n 
asit’};  the  .Vtcrosoft  Opt  i.m  i /.  i  ny  FORTRAN  compiler  version  4.10.  The 

fiutales  wt;re  eominleii  asinji  The  Microsoft  QuickC  co!;’[)iler  version  1.00. 
fin;  ASSr.MHl.i'.'R  modal (.‘s  wer<;  assembled  with  the  Microsoft  .Macro  Asst.’wbler 
ve.’-sion  ,■).().  Tht;  motiules  were  linked  with  the  Microsoft  Extended  Overlay 
i.tttker  version  5.01.20  ana  all  software  was  developed  on  an  TB.M  PC/XT 
cmii’paiiiin;  atnu;r  the  PC-I)0S  version  3.1  op(;ratinij  system. 


Qi. 
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Si  final  Proci.-ssiiifi  for  HC!  Monitor. 
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Fiftare  B.l. 


Figure  B.2.  FJowcliart  for  HCl  Program 


SOL'RL'i;  CODE  FOR  llCl  FROGRiVM 


c 

c...  Suoroutine  to  analyze  the  voltage  array  A  from  the  main 
c...  HCl  program.  The  major  result  is  a  returned  value  of  X 
c...  which  is  related  to  the  HCl  concentration  in  ppm.  If  full 
c...  output  is  requested,  as  in  a  calibration  run,  then  a  number 
c...  intermediate  results  are  displayed  to  the  screen, 
c 

subroutine  analOl (f lag,n,rO,Rtot ,X,A) 
c 

c...  A(i,j)  =  Array  containing  the  voltage  measurements  of  all 

c. . .  of  the  wheel  turns 

c...  Ared(i)  =  peak-null  voltage  of  first  peak 

c...  Bred(i)  =  peak-null  voltage  of  second  peak 

c...  Cred(i)  =  peak-null  voltage  of  third  peak 

c...  Ored(i)  -  peak-null  voltage  of  fourth  peak 

c...  Oiffl(i)=  {1st  peak  -  2nd  peak  difference)/8red 

c...  0iff2(i)=  (3rd  peak  -  4th  peak  uifference)/Ored 

c...  R(i)  =  HCl  ratio  for  a  single  wheel  turn 

c...  sigma(i)=  standard  deviation  of  ave(i) 

c...  ave(i)  =  average  of  the  eight  measurements,  Ared,  Bred,  Cred,  Ored 
c...  Oitfl,  and  Diff2  over  the  number  of  wheel  turns 

c. . .  Rtot  =  HCl  ratio 

c...  X  =  -ln(Ratio)  related  to  HCl  concentration 

c...  n  =  number  of  wheel  turns 

c...  rO  =  HCl  ratio  in  the  absence  of  HCl 

c 

real*4  A(8,401),Ared(400),Bred(400),Cred(400) 
real*4  Dred(400),Diff 1 (400), Oiff 2(400), R(400) 
real*4  sigma{14),ave(14),Rtot,X 
character  flag’l 
c 

it{ .not. (tlag.eq. 'y' .or.flag.eq. 'Y' .or.flag.eg. 'n' 

1  .or.flag.eq. 'N' ))  then 
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print  ♦.’Input  error  to  analOl' 
print  *, 'Flag  must  equal  "Y*  or  ‘N’.  input  flag' 
stop 
end  if 

Rtot  -  0.0 
do  10  j=l ,n 

Ared(j)  -  A(2,j)  -  0.5»(A(1,j)  +  A{3,j)) 

8red(j)  =  A(4,j)  -  0.5*(A(3,j)  +  A(5.j)) 

Cred(j)  =  A(6,j)  -  0.5»(A(5,j)  +  A(7,j)) 

Ored(j)  =  A(8,j)  -  0.5»(A(7.j)  +  A(1,(j+1))) 

Diffl(j)  -  (A(2.j)  -  A(4,j))/3red(j) 

Ditf2(j)  =  (A(6,j)  -  A(8.j))/0red(j) 

R(j)  =  (A(2,j)*A(6.j)-A(4,j)-A(8,j))/(Bred(j)i0red(j)) 

Rtot  =  Rtot  +R(j) 

continue 

Rtot  =  Rtat/(r0»f loat(n)) 
i f (Rtot. le. 0.0)  Chen 
X  =  5.0 
else 

X  =  -alog(Rtot) 
end  if 

do  100  i=l,8 

ave(i)=0.0 
do  50  j=l ,n 

ave(i)=ave(i)+A(i,j) 

continue 

continue 

ave(9)=0.0 
ave(10)=0.0 
ave(11)=0.0 
ave{12)=0.0 
ave( 13)=0.0 
ave{ 14)=0.0 


do  110  j=1,n 


110 

c 

c 


150 

200 

c 


ave(9)  =  dve(9)  +  Ared(j) 
ave(IO)  =  ave(IO)  +  8ped(j) 

ave(11)  =  dve(11)  +  Cred(j) 

ave(12)  =  ave(12)  +  Dred(j) 

ave(13)  =  ave(13}  +  Oiffl(J) 

ave(14)  =  ave{14)  +  Diff2(j) 


continue 


do  200  1=1,8 

ave(1)=ave(1)/float(n) 
sl9ma(1)=0.0 
do  150  j=1,n 

s1gma(1)=s1gma(1)+(A(1,J)-ave(1))**2 

continue 

s1gnia(1)=sqnt(s1gtna(1)/f1oat(n)) 

continue 

ave(9)=ave(9)/float(n) 

ave(10)=ave(10)/float(n) 

ave(11)=ave(11)/f1oat(n) 

ave(12)=ave(12)/float(n) 

ave(13)=ave(13)/float(n) 

ave(14)=ave(14)/tloat(n) 

s1gma(9)=0.0 

s1gma(10)=0.0 

s1gina(11)=0.0 

s1gma(12)=0.0 

s1gma(13)=0.0 

s1gnia(14)=0.0 

do  250  j=1  .n- 

s1gma(9)=s1gma(9)+(Ared(J)-ave(9))**2 
s1gma(10)=s1gma(10)+{Bred(j)-ave(10))**2 
sigma(  1 1  )=s1giiia(  1 1  )+(Cred(  j)-ave(1 1  ))»*2 
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sigma(12)=sigitia(12)+(0red(J)-ave(12))**2 
sigma(13)=sigma(13)+(Diff l(J)-ave(13))**2 
sigma(14)=sigiBa(14)+(D1ff2(J)-ave(14))**2 

250  continue 

c 

do  270  1=9,14 

s1gina(1)=sqpt(s1gma(1)/float(n)) 

270  continue 

c 

If (f lag.eq. 'Y' .or.f1ag.eq. 'y')  then 
c 

wr1te(6,*)'  ' 

»/r1te(6,*)'  ' 

wr1te(6, *)'****»*  Statistics  concerning  raw  data  »»»»»»’ 

wr1te(6,*)'  ' 

wr1te(6,ni) 

111  foriiiat(11x,'r.8x.'2',8x,'3’,8x,'4',8x,’5',8x/6’,8x,'7',8x,'8') 
wr1te(6,1)  (ave(1),1=1,8) 

1  forinat(1x,'ave',3x,F8.5.7(1X,F8.5)) 
wr1te(6,2)  (s1giiia(1),1=1,8) 

2  format(1x,'s1gma',1x,F8.5,7(1X,F8.5)) 
wr1te(6,*)'  ' 

wr1te(6,*)'  ' 
wr1te(6,*)'  ' 
wr1te(6,44) 

44  format( lx, '*♦♦♦♦*  Statistics  concerning  four  peaks', 

r  and  difference  ratios  »»**»»') 
wr1te(6,*)'  ' 
wr1te(6,112) 

112  f ormat { 1 X , 9x , ' Peaki ' , 7x , ' Peak2 ' , 7x , ' Peak3 ' , 7x , ' Peak4 ' , 

1  7x,'D1ff1',7x,'D1ff2') 

wr1te(6,113)  (ave(1),1=9,14) 

113  format (lx, 'ave' ,3x,F8.5,5(4x,F8.5)) 
wr1te(6,3)  (s1gma(1),1=9,14) 

3  format (lx, 'sigma' , 1x,F8.5,5(4x,F8.5)) 
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c 


wpite{6,*)'  ' 

wn'te(6,»)'The  swaller  difference  ratio  average  1s', 


1  atninl(ave(13),ave(U)) 

write(6,*)'The  larger  difference  ratio  average  1s', 
1  amaxl(ave(13),ave(U)) 
endif 

tet!!l=(ave(  13)+ave(U))/2.0 
if(flag.eq.'Y'.or.flag.eq.'y')  then 
write(6,*)'The  average  of  averages  is',teffl1 
endif 

teml=aminl (ave(13) ,ave(14)) 

terr,2=amaxl{ave(13),ave(14)) 

tem3=(ave(13)+ave( 14))/2.0 

return 

end 

interface  to  subroutine  gdate  (x) 
inteqer’4  x 
end 
c 

subroutine  getdat (month, day, year) 
c 

c  Routine  to  get  the  date  through  an  MS-DOS  int21H 
c  function  call 
c 

integer’4  year, month, day, date 
c 

call  gdate(date) 
c 

year  =  iand(date, 16#FFFFOO00)/16#1O0O0 
month  =  iand(date, 16#FFOO)/16#lOO 
day  =  iand(date, 16#FF) 
c 

return 

end 
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interface  to  subroutine  gtime  (x) 
integer*4  x 
end 
c 

subroutine  gBttim(hours, minute, second, hundre) 
c 

c  Routine  to  get  the  time  of  day  through  an  MS-DOS  int21H 
c  function  call 
c 

integer*4  hours, minute, second, hundre, time 
c 

call  gtime(time) 
c 

hours  =  iand(time,  16tIFFOO0000)/16tl1O00OOO 
minute  =  iand(time,  16#FFQQOO)/16ltlOOOO 
second  =  iand(time, 16#FF0O)/16#1OO 
hundre  =  iand(time, 16#FF) 
c 

c  print  101 , hours, minute, second, hundre 

101  format(lx,i2,':',i2,':',i2,'.',i2) 

return 
end 

interface  to  subroutine  pltgr  [c]  (x) 

integer*4  x  [near,  reference] 

end 

interface  to  subroutine  pltpts  [c]  (x,y,z) 
real*4  x  [near,  reference] 
real*4  y  [near,  reference] 

1nteger*4  z  [near,  reference] 
end 

interface  to  subroutine  rscrn  [c] 
end 
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interface  to  subroutine  eschit  [c]  (x) 

integer*^  x  [near,  reference] 

end 


c 

c...  This  is  the  main  program  unit  for  the  SSI  HCl  mor-itor 
c 


program  hcl 
c 

real*4  A(8,401 ) .tinit ,tf inal ,etim,cppm 
real*4  rO,rOa,T 

integer*4  Yval(3208),numave,ngain, lent, mtype, heat 
integer»4  hrs , min , sec , hun , hkey ,ncal 
character  flag^l ,filnam*'14 
c 


c. . 

•  A(i,j) 

=  Array  in  which  the  measured  voltages  are  stored 

c. . 

.  Yval(i) 

=  Array  which  contains  the  integers  of  the  A/D  conversions 

c. . 

.  tinit 

=  Initial  time 

c. . 

.  tfinal 

=  Final  time 

c. . 

.  etim 

=  Elapsed  time 

c. . 

.  cppm 

=  HCl  concentration  in  ppm 

c. . 

.  rO 

=  Measured  ratio  in  the  absence  of  HCl 

c. . 

.  rOa 

=  Dummy  variable  used  in  calibration 

c. . 

.  T 

=  Ambient  temperature  in  Celsius 

c. . 

.  numave 

=  Number  of  wheel  turns  to  average 

c. . 

.  ngain 

=  Gain  for  which  the  A/0  board  is  set 

c. . 

.  lent 

=  Length  of  time  for  measurement  in  minutes 

c. . 

.  mtype 

=  Measurement  type,  calibration(O)  or  HCl  determination(l) 

c. . 

.  heat 

=  Aerosol  heater  on(1)  or  off(O) 

c. . 

.  hrs 

=  Hours 

c. . 

.  min 

=  ''inutes 

. 

.  sec 

=  Seconds 

c. . 

.  hun 

=  Hundredths  of  seconds 

c. . 

.  hkey 

=  Value  is  27  is  escape  key  hit,  0  otherwise 
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c...  ncal  =  Increment  tor  calibrtion  measurement 
c...  filnam  =  File  name  in  which  to  store  measurements 
c 
c 

c...  cone  is  the  function  which  relates  the  In  of  the  measured 
c...  ratio  and  the  ambient  temperature  to  the  HC1  concentration 
c. . .  in  ppm 
c 

conc(T,X)  =  ({273.0+T)/295.0)»*0.75»{95.08*X  ♦  20.59»X»X  + 

&  19.65»x»»3  -  1.893»X*»4  +  0.3074»X»*5  +  0.2290»X**6) 
c 

c...  hclscr  sets  up  the  initial  screen  and  holds  it  for  3  seconds 
c 

call  hclscr(3.0) 
c 

c...  rdparm  is  the  subroutine  which  reads  all  of  the  user  input  and 
c...  relays  the  information  back  to  the  main  program 
c 

call  rdparm(numave, ngain.f actor, lent,filnam,rO,mtype,heat,T) 
c 

c...  If  mtype=1  then  do  an  HCl  measurement 
c 

if(mtype.eq.l)  then 
c 
c 
c 

call  gettim(hrs,min,sec,hun) 
tinit  =  3600*hrs  +  60*min  +  sec  +  hun/100 
c 

c...  pltgr  is  the  subroutine  which  sets  up  the  graphics  screen 
c 

call  pltgr(lent) 
c 

c...  stpclr,  setad,  and  rdAtoD  set  up  and  read  the  A/0  board 
c 
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10  call  stpclr 

ca 1 1  setad(ngain , 8*nufflave+7 ,0,0) 
call  rdAtoO(8*numave+7, 'e', 'i',Yval) 
c 

c...  sorter  sorts  Yval  according  to  tne  proceourt  outlined  in  appendix  8 
c...  of  the  final  report 
c 

call  sorter(nuroave,Yval) 
c 

c...  convert  Yval  to  volts  array  A 
c 

do  100  j  =  1,nuinave+1 
c 

do  1  =  1,8 

A(iJ)  =  ((20.0  *  Yva1(b*(j-l)+i)/4096.0)  -  10.0)/factor 

44  continue 

c 

100  continue 

c 

c...  analyze  the  array  A  for  the  In  of  the  average  ratio  X 
c 

call  anal01('n' ,numave,r0,Rtot,X,A) 
call  gettin(hrs,niin,sec,hun) 
tfinal  =  3600*hrs  +  SO^ntin  +  sec  +  hun/100 
etim  =  (tf inal-tinit)/60.0 
if(X.ge.5.0)  then 
cppin  =  10000.0 
else 

cppm  =  conc(T,X) 
end  if 
c 

c...  plot  the  HCl  concentration  and  store  in  filnara 
c 

call  pltpts(etim, cppm, lent) 
write(20,*)  etim, cppm 
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c 

c...  check  if  measurement  time  is  completed  or  the  escape  key  is  hit 
c...  if  yes  close  filnam  and  clear  and  reset  screen  to  text  mode 
c...  if  no  continue  taking  data 

r 

c 

if(etim.ge.f loat(lent))  then 
close(unit=20) 
call  rscrn 
stop 
else 

call  aschitfhkey) 
if(hkey.eq.27)  then 
close(unit=20) 
call  rscrn 
stop 
else 
goto  1C 
end  if 
end  if 
c 
c 
c 

end  if 
c 
c 

c...  If  mtype=0  then  do  a  calibration  measurement 
c 

if (mtype.eq.O)  then 
c 
c 
c 

rOa  =  0.0 
rO  =  t.O 
ncal=l 
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print  *, 'Calibration  step' 
print  ' 
c 

call  gettiiii(hrs,min,sec,hun) 
cinit  ^  3600»h.o  +  oO*iittn  +  sec  +  hun/100 
c 

c...  stpcir,  setad,  and  rdAtoO  set  up  and  read  the  A/0  board 
c 

20  call  stpcir 

call  setafl(ngain,8*nuniave+7,0,0) 
call  rdAtoO(8»numave+7.'e','i'.Yval) 
c 

c...  sorter  sorts  Yval  according  to  the  procedure  outlined  in  appendix  B 
c. . .  of  the  final  report 
c 

call  sorter(numave,Yval) 
c 

c...  convert  Yval  to  volts  array  A 
c 

do  200  j  =  l,numave+1 
c 

do  24  i  =  1,8 

A(1,j)  =  ((20.0  »  Yval(8»(j-l)+i1/4096.0)  -  10.0)/factor 
24  continue 

c  write(25,234)  (A(i,j),i=1,8) 

234  forinat(lx,8t9.4) 

c 

200  continue 

c 

c...  analyze  the  array  A  for  the  In  of  the  average  ratio  X 
c 

call  anal01( 'y' ,numave,rO,Rtot,X,A) 
call  gettim(hrs,min,sec,hun) 
tfinal  =  3600»hrs  +  60*inin  +  sec  +  hun/100 
etim  =  (tfinal-tinit)/60.0 
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cppm  =  conc(T,X) 


c 

c. . .  store  X  In  file  filnam 
c 

write(20,*)  etira.X 
c 
c 

c...  take  punnig  average  of  new  rO 
c 

rOa  =  rOa  +  exp(-X) 
ncal=ncal+l 
c 

c...  check  if  measurement  time  is  completed  or  the  escape  key  is  hit 

c...  if  yes  close  filnam,  store  new  rO  in  file  calib, 

c...  and  clear  and  reset  screen  to  text  mode 

c...  if  no  continue  taking  data 

c 

c 

if (etim. ge. float (lent ). or. ncal.gt. 500)  then 
close(unit=20) 
rOa  =  r0a/float(ncal-1) 
rewind(unit=15) 
write{15,*)  rOa 
close(unit=15) 
call  rscrn 
stop 
else 

call  eschit(hkey) 
if (hkey.eq.27)  then 
close(unit=20) 
rOa  =  rOa/float(ncal-l) 
rewind(unit=15) 
write(15.*)  rOa 
close(unit=15) 
call  rscrn 
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stop 
else 
goto  20 
end  if 
end  it 
c 
c 
c 

endif 

c 

c 

c 

1000  stop 
end 

interface  to  subroutine  els 
end 

interface  to  subroutine  pstrng  (x) 

character  x»81 

end 

subroutine  hclscr(second) 
c 

c  ...  This  subroutine  clears  the  screen  and  prints  the  SSI 
c  ...  copyright  information  and  leaves  it  up  for  the  number 
c  ...  of  seconds  specified.  Then  the  screen  is  cleared 
c  . . .  once  more, 
c 

integer*4  i 
real ’*4  second 

character  c*81 ,str1»40,str2*40 
c 

call  els 

c  write(6, '(lx,80(lH*))') 

c(  1 :40)=str1 
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c(41  :80)=StPl 
C(81:81)='$' 
call  pstrng{c) 
str1='» 

str2=’  *■ 

c(  1 :40)=strl 

c(4).-80)=str2 

c(81:81)='$' 

do  10  i = 1 , 3 

c  write(6,'(1x,1H*,78(lH 

call  pstrng(c) 

10  continue 

strl='»  HCl  MO’ 

str2='NITOR  »- 

c( 1 ;40)=str1 

c(41;80)=str2 

c(8l:81)='{' 

call  pstrng(c) 

c  write(6, '(lx,1H»,34(lH  ),11HHC1  MONITOR, 33(1H 

strl='*  • 

str2=-'  »' 

c(  1 :40)=str1 
c(41;80)=str2 
c(81:81)='$' 
do  20  1=1,3 

c  wpite(6,'(1x,lH»,78(1H  ),1H*)‘) 

call  pstPng(c) 

20  continue 

c  wpite(6, '  (lx,  1H*,  19(lil  ) ,40Hcopypight  1938,  Spectpal  Sciences,  Inc., 
c  &  19(1H  ),1H*)') 

stPl='*  copyright  1988,  Spe' 

stp2='ctpal  Sciences,  Inc.  *' 

c(1  •.40)=stp1 
c(41 :80)=stp2 
c(81;81)='$' 
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call  pstpng(c) 
strl='* 

str2='  *' 

c(1:40)=str1 

c(41  ;80)=str2 

c(81;81)='$' 

do  30  1=1,9 

c  wnte(6,’(1x,1H«,78(lH 

call  pstrng(c) 

30  continue 

c  write(6,'(1x.80(1H»))') 

str 1= ' **************************************** ' 

c(1;40)=str1 

c(41:80)=str1 

c(81:81)='r 

call  pstrng(c) 

call  ttta It (second) 

call  els 

return 

end 

interface  to  subroutine  inplb  (a,b) 

integer*2  a,b 

end 

interface  to  subroutine  outlb  (a,b) 

integep*2  a,b 

end 

SUBROUTINE  RdAto0(nreads,ieclk,1etr1g,Yval) 
c 

c... FORTRAN  routine  to  read  the  OT-2401  board  following  A/D  setup  with 
c,. .‘setup*.  The  inputs  are:  number  of  conversions,  flag(i,e)  for  internal 
c...or  external  clock,  and  flag(i,e)  for  interna!  or  external  trigger. 
c...The  output  is  an  array  of  A/0  conversions 
c 
C 
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c 

C  6SA0RS  =  Base  Address 

C  CMOREG  =  Command  Register 

C  CMDWAI  =  Command  Wait 

C  DATREG  =  Data  Register 

C  ROWAI  =  Read  Wait 

C  STAREG  =  Status  Register 

C  WRTWAI  =  Write  Wait 

C  AOL  =  Array  to  hold  low  Byte  of  data 

C  ADH  =  Array  to  hold  high  Byte  of  data 

C  EXTCLK  =  External  Clock  code 

C  EXTRIG  =  External  Trigger  code 

C 

integer*^  nreads 
integer*4  Yval(3208) 
character  iec1k*1 ,ietrig*l 
C 

I NTEGER*2  BSAORS , CMDREG , STAREG , DATREG , CMOWAl , WRTWAI 
INT£GER»2  ROWAI .CSTOP, STATUS, EXTRIG. ERR0R1 ,ERR0R2 
INTEGER»2  CRAO,ADL(4008),AOH(4008), EXTCLK, COMAND 
C 
C 

10  BSADRS  =  748 

CMDREG  =  BSAORS  +  1 
STAREG  =  BSADRS  +  1 
DATREG  =  BSADRS 
CMDWAI  =  4 
WRTWAI  =  2 
RDWAI  =  5  . 

C 

C 

CSTOP  =  15' 

CRAD  =  14 

EXTCLK  =  64 
EXTRIG  =  128 
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c 

c 

C  Check  for  legal  Status  Register  value. 

C 

CALL  INPLB{STAREG, STATUS) 
c 
c 

IF  (.NOT.((IAND(STATUS,112)).EQ.O))  60  TO  98 
C 
C 

COMANO  =  0 

i f ( iec 1 k . EQ . ' E ' . OR . 1 ec 1 k . EQ . ' e ' )  COMANO=COHAND+EXTCLK 
if{ietrig.E0.'E'.OR.letrig.EQ.'e')  C0MAN0=C0HAND+EXTR1G 
C 

comand  =  crad  *  comand 

call  HAIT(STARE6,CHOWA1.0) 
call  OUTLB(CMDREG, comand) 

C 

DO  222  LOOP  =  l.nreads 

CALL  WA1T(STAR£6,R0WAI,0) 

CALL  INPL8{DATREG,A0L(L00P)) 

CALL  «AIT(STAREG,R0WAI,0) 

CALL  INPLB(OATREG,AOH(LOOP)) 

222  CONTINUE 

C 

C  Check  for  ERROR. 

C 

CALL  WA1T(STAREG,CMOWA1,0) 

CALL  INPLB(STAREG, STATUS) 

C 

c  I  added  this  line  in  order  to  ignore  the  high  4  bits 

c 

c 

STATUS  =  IAN0(STATUS,15) 
c 
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IF(1AN[)(STATUS,128).NE.0)  GOTO  94 


C 

DO  333  LOOP  =  l.nreads 
YvaULOOP)  =  ADH{LOOP)*256  +  ADL(LOOP) 
if{Yva1(L00P).6T. 32767)  Yva](LOOP)  =  Yval(LOOP)  -  65536 
333  CONTINUE 

C 

69  continue 

c  WR1TE(*,»)'  ' 

c  WR1TE(*,»)'  ' 

c  WRITE(*,*)  '  READ  A/O  Operation  Complete’ 

GOTO  102 
C 

C  Fatal  board  error. 

C 

94  WRITE(*,»)  '  ' 

HR1TE(*,»)  'FATAL  BOARD  ERROR' 

WR1TE(6,2)  STATUS 

2  F0RMAT( IX, 'STATUS  REGISTER  VALUE  IS  ',12,'  DECIMAL’) 

MRITE(»,«)  '  ' 

CALL  SR9600 (CM0RE6 , STAREG , DATREG , CMOWAI . ROWAI . 


-CST0P,CERR0R,ERR0R1,ERR0R2) 

WRITE!*,*)  'ERROR  REGISTER  Vt:,JES  ARE:' 
WRITE(6,3)  ERRORl 


3 

FORMAT! IX,'  BYTE  1  - 

',12,' 

DECIMAL') 

WRITE!6,4)  ERROR2 

4 

FORMAT! IX,'  BYTE  2  - 

',12,' 

DECIMAL’) 

WRITE!*,*)  '  ' 

GOTO  102 
C 

C  Illegal  Status  Register. 

C 

98  WRITE!*,*)  '  ' 

WRITE!*,*)  'FATAL  ERROR  -  ILLEGAL  STATUS  REGISTER  VALUE' 
WR1TE!6,5)  STATUS 
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5 


FORMAT( IX, 'STATUS  REGISTER  VALUE  IS  '.12,'  DECIMAL') 
write(*,*)'tailure  in  ‘rdAtoO"' 

C 

102  RETURN 
END 

interface  to  subroutine  inplb  (a.b) 

integer*4  a.b 

end 

interface  to  subroutine  out'lb  (a.b) 

integer*4  a.b 

end 

subroutine  setad(GNCOOE,NUMCNV,SRTCHN.ENDCHN) 

C 

integer*4  STAREG.STATUS.CMDWAI .CMOREG.CSAO.WRTMAI 
integer»4  DAT  REG , 6NC0DE , SRTCHN , ENDCHN . BSADRS 
integer^  NUML.NUHH.NUMCNV 
c 

c... FORTRAN  routine  to  set  the  DT-2401  board  for  the  proper  A/D  gain, 
c... number  of  conversions,  start  channel,  and  end  channel.  To  be  called 
c... following  'stpclr* 
c 
C 
C 

C  SSADRS  =  Base  Address 

C  CMDREG  =  Command  Register 

C  CMDWAI  =  Command  Wait 

C  OATREG  =  Data  Register 

C  STAREG  =  Status  Register 

C  KRTWAI  =  Write  Wait 

C  GNCOOE  =  Gain  Code 

C 
C 

10  BSADRS  =  748 

CMDREG  =  SSADRS  +  1 


52 


STAREG  =  BSADRS  +  1 
DATREG  =  BSADRS 
CMDWAI  «  4 
HRTWAl  =  2 
C 
C 

CSAD  =  13 
C 
C 

C  Check  for  legal  Status  Register  value. 

C 

CALL  INPLB(STAREG, STATUS) 
c 
c 

IF  (.NOT.((IANO(STATUS,112)).EQ.O))  GO  TO  98 
C 

C  Do  a  SET  A/D  PARAMETERS  command  to  set  up  the  A/D  converter. 
C  Write  SET  A/0  PARAMETERS  command. 

C 

CALL  WAIT(STAREG. CMDWAI, 0) 

CALL  OUTLB(CMOREG,CSAO) 

C 

C  Write  A/O  gain  byte. 

C 

CALL  WAIT(STAREG,WRTWAI,WRTWAI) 

CALL  OUTLB(DATREG,GNCODE) 

C 

C  Write  A/0  start  channel  byte. 

C 

CALL  WA1T(STAREG,WRTWA1,WRTWAI) 

CALL  OUTLB(DATREG,SRTCHN) 

C 

C  Write  A/D  end  channel  byte. 

C 

CALL  WAIT(STAREG,WRTWA1,WRTWA1) 
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CALL  OUTLB(OATREG,ENDCHN) 


C 

C  Write  high  and  low  bytes  of  NCONVERSIONS#. 

C 

NUMH  =  lNT(NUMCNV/256) 

NUML  =  NUMCNV  -  NUMH  »  256 
CALL  HAIT{STARE6,WRTWA1,HRTWAI) 

CALL  OUTL8(DATREG,NUML) 

CALL  WAIT(STAREG,WRTWAI.WRTWAI) 

CALL  OUTLB(OATREG,NUMH) 

C 

c  write(*,*)'  ' 

c  write(».»)'SET  A/D  PARAMETERS  completed’ 

goto  102 

98  WRITE(*.»)  '  ' 

WRITE(*,*)  'FATAL  ERROR  -  ILLEGAL  STATUS  REGISTER  VALUE’ 
HRITE(6,5)  STATUS 

5  FORMAT( IX. ’STATUS  REGISTER  VALUE  IS  ’.12,’  DECIMAL') 

write(*,*) ’failure  in  'setad" 

C 

102  return 

end 

subroutine  sorter(n,Y) 
c 

c... FORTRAN  routine  to  sort  the  data  read  from  the  DT-2401  board  so 
c...that  the  first  point  is  a  baseline,  the  second  is  the  first  large 
c... maximum,  the  third  is  a  baseline,  the  fourth  is  the  first  small 
c... maximum,  then  a  repetition  for  up  to  the  number  of  wheel  turns 
c. . .specified, 
c 
c 

integer*4  n,'Y(3208) , dum, ndum, idum.mdum 
c 


dum  =  max(Y(l),Y(2),Y{3),Y(4),Y(5),Y(6),Y(7)) 


10 

c 

11 

12 

c 


13 
c 

14 


20 


c 


30 


mdum  =  -32767 

do  10  i=1,7 

1t(Y(i).EQ.dum)  then 
ndum  =  i 
goto  11 

end  if 

continue 

do  12  i=1,7 

if (i .NE.ndum)  mdum  =  max(mdum,Y(i)) 

continue 

do  13  i=1,7 

if (Y(i) .EO.mdum.AND. i .NE.ndum)  then 
kdum  =  i 
goto  14 

end  if 

continue 

if (ndum. GE. 2)  then 

do  20  i=l,8»n+1 

Y(i)  =  Y(i+nduro-2) 

continue 

return 


else 

if (kdum. GE. 2)  then 

do  30  i=1,8*n+1 

Y(i)  =  Y(i+kdum-2) 

continue 

return 


else 

write(*,*) error  in  subroutine  sorter 
stop 
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end  if 


end  it 

return 

end 

interface  to  subroutine  inplb  (a,b) 

inteqer*4  a,b 

end 

interface  to  subroutine  outib  (a.b) 

integer*4  a,b 

end 


SUBROUTINE  SR9600(CMDRE6,STAR£G,OATR£G.C«OHAI,RDWAI, 
-CST0P,CERR0R,ERR0R1,ERR0R2) 

c  This  part  should  be  separated  into  a  separate  subroutine 
c 

c... FORTRAN  routine  to  read  messages  from  the  DT-2401  board 
c 

INTEG£R*4  CMDREG . STAREG , OATREG .CMOWAI , ROWAl 
INTEGERS  CST0P,C£RR0R,TEMP,ERR0R1,ERR0R2 
C  Read  the  Error  Register. 

C 

CALL  OUTLB (CMDREG, CSTOP) 

CALL  INPLB(DATREG,TEMP) 

C 

CALL  HA1T(STAREG,CMDWAI,0) 

CALL  OUTLB(CMDREG,CERROR) 

C 

CALL  WA1T{STAREG,RDHAI,0) 

CALL  INPL8{DATREG,ERR0R1) 

C 

CALL  WA1T(STAREG, ROWAl, 0) 

CALL  INPLB(DATREG,ERR0R2) 

C 

RETURN 
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END 

interface  to  subroutine  inplb  (a,b) 

integer*4  a,b 

end 

interface  to  subroutine  outlb  (a,b) 

integerM  a,b 

end 

SUBROUTINE  stpc'lr 
c 

c... FORTRAN  routine  to  stop  and  clear  the  DT-2401  board 
c 
C 
C 

C  6SADRS  =  Base  Address 

C  CMOREG  =  Command  Register 

C  CMDWAI  =  Command  Wait 

C  0ATRE6  =  Data  Register 

C  STAREG  =  Status  Register 

C 

c 

INTEGER*4  BSADRS, CMOREG, STAREG, OATREG.CMOHAI 
INTEGER»4  TEMP, CCLEAR,CSTOP, STATUS 
C 

BSADRS  =  748 

CMOREG  =  BSADRS  +  1 

STAREG  =  BSADRS  +  1 

OATREG  =  BSADRS 

CMDWAI  =  4  . 

C 

C 

CCLEAR  =  1  ■ 

CSTOP  =  15 

C 

C  Check  for  legal  Status  Register. 
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c 

CALL  INPL8(STARE6. STATUS) 
c 
c 

1F(.NOT.((IANO(STATUS,112)).EQ.O))  GOTO  98 
C 
C 

C  Stop  and  clear  DT2801. 

C 

CALL  OUTLB(CMDR£G.CSTOP) 

CALL  INPLB(OATREG,TEMP) 

C 

CALL  WAIT(STAREG,CMDWA1,0) 

CALL  OUTLB(CHOREG,CCLEAR) 

C 

c  write(*,*)'STOP  and  CLEAR  completed' 

goto  102 
C 

C  Illegal  Status  Register. 

C 

98  WRIT£(*,»)'  ' 

HRITE(»,*)  'FATAL  ERROR  -  ILLEGAL  STATUS  REGISTER  VALUE' 
WRITE(6.2)  STATUS 

2  FORMAT! IX, 'THE  VALUE  OF  STATUS  IS  ',12,'  DECIMAL') 

C 

102  RETURN 

END 

subroutine  twait(s) 
c 

c... FORTRAN  routine  to  wait  a  certain  amount  of  time  (sec)  then  precede 

c 

c 

real*4  s,tinit,tf inal 
integer*!  i , n , i hrs , imin , isec , i hun 
c 
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call  gettim(ihrs,imin,isec,ihun) 
tinit  =  f1oat(3600*ihrs  +  60*imin  +  isec  +  0,01*ihun) 
call  gettim( ihrs , Imin , isec , ihun) 
tfinal  =  f1oat(3600*ihrs  +  60*imin  +  isec  +  0.01*ihun) 
if((tfinal-tinic).lt.s)  goto  10 
return 
end 

interface  to  subroutine  inplb  (a,b) 
integer*4  a,b 
end 

subroutine  wait(addres,n,m) 
c 

c... FORTRAN  routine  to  suspend  execution  until  a  specified  bit  pattern 
c...is  read  from  the  specified  port 
c 

integer*4  addris, value 
integer*4  result, n,m 
1  call  inplb(addres, value) 
result  =  I£OR(value,m) 
result  =  ;ANi:;result,n) 
if(result.eq.O)  go  to  1 
return 
end 
c 

c...  Subroutine  to  read  all  of  the  user  supplied  parameters  for  the  HCl 
c...  program, 
c 

subroutine  rdparm(numave, ngain.f actor, lent,filnam,rO,mtype, heat ,T) 
c 


c. . 

.  nunave 

=  Number  of  wheel  turns 

to  measure  over 

c. . 

.  nqain 

=  Gain  of  A/0  board 

c. 

=  Oivisicn  factor  based 

on  ngain  tor  voltage  conversion 

c, . 

.  lent 

=  Time  length  for  measurement  run 

c. . 

.  filnam 

=  Name  of  user  supplied 

file  to  store  HCl  concentrations 
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c...  rO  =  HCl  ratio  in  absence  of  HCl  absorption 

c...  mtype  =  Measurement  type;  0  for  calibration,  1  for  HCl  determination 

c...  heat  =  Aerosol  heater  on(1)  or  off(O) 

c...  T  =  Ambient  tempprature  in  Celsius 

c 

integer’A  numave.ngain, lent , mtype, heat 
real*4  rO, factor,! 
character  filnam*14,header’'72 
c 

write(6,*) 'Average  over  how  many  measurements  ' 
read(5,*)  numave 

c33  write(6,*) 'select  gain  factor:' 
c  write(6,*)  '(0  for  gifl  V)' 

c  write(6.’)  '(1  for  g5  V)' 

c  write(6,*)  '(2  for  g2.5  V)' 

c  write(6,*)  '(3  for  gl .25  V)' 

c  write(6,*) 'gain  factor:  ’ 

c  ''®ad(5,’^)  ngain 

c  if (ngain.LT.O.OR.ngain.GT.B)  then 

c  write(6,’*)'  ' 

c  write(6,*) 'Please  input  a  number  between  0  and  3' 

c  goto  33 

c  end  if 

ngain=2 
write(6,*)'  ' 
if(ngain.eq.O)  factor=1.0 
if(ngain.eq.l)  factor=2.0 
if (ngain. eg. 2)  factor=4.0 
if (ngain. eg. 3)  factor=8.0 
c 

5  write(6,*)'Name  the  file  you  wish  to  store  data  to' 

11  format(Al) 

12  format(A14) 

13  format(A72) 

14  fornat( lx,A72) 
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c 

write{6,’')'tilenaiiie=' 
i'ead(5,1?)  filnatn 

write(6,*)'print  a  one  line  header' 
reaa(5,13)  header 

cpen(20,f  ile=tilnaiii,  status^' unknown ' ) 
write(20,14)  header 

0 

mtype  =  -1 

50  write(6,*)  'Measurement  type' 

write(6,*)  'Calibration  =  O' 
write(6,*)  'HCl  testing  =  1 ' 
read(5,*)  mtype 

if(. not. (mtype. eq.O. or. mtype. eg. 1))  goto  50 
heat=-; 

55  write(6,*)  'Aerosol  measurement' 

write(6,*)  'Aerosol  Heater  off  =  O' 
write(6,'')  'Aerosol  Heater  on  =1' 
readfS,*)  heat 

write(6,*)  'Input  the  ambient  temperature  in  Celsius' 
read(5,*)  T 

if(. not. (heat. eq.O. or. heat. eg. 1))  goto  55 
write(6,*)  'Input  time  duration  of  data  accumulation' 
readlS,*)  lent 

open(unit=15,f ile='calib' ,status='old' ) 

if (mtype. eg . 1 )  read (15,*)  rO 

return 

end 

llinclude  <stdio.h> 
ttinclude  <conio.h> 
tinclude  <ctype.h> 

/*  This  C  function  checks  whether  or  not  the  escape  key  was  hit 
yes  then  27  (ascii  representation  of  escape  key)  is  returned 
no  0  is  returned  */ 
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/*  called  as  FORTRAN  subroutine  ‘eschitCdum)' 
integer*4  dum 

void  eschit(duffl) 
int  near  ♦dum; 


int  c; 

if  (  kbhitO  )  ( 
c  =  getchO: 

} 

else  { 
c  =  0; 

1 

if  (c==toascii(27))  ♦dum  =  27; 
if  (c!=toascii(27))  ♦dum  =  0; 

) 

Sinclude  <stdio.h> 

((include  <ctype.h> 

((include  <graph.h> 

struct  videoconfig  vc;  /♦  variable  vc  of  type  videoconfig  ♦/ 
void  pltgr(tmax) 

int  near  ♦tmax;  /*  tmax  is  the  time  scale  (min)  maximum  ♦/ 

/♦  This  C  function  sets  up  the  plotting  coordinates  for  the  HCl 

real  time  measurements.  Used  with  'pltpts'  to  plot  the  points  ♦/ 
/♦  called  as  FORTRAN  subroutine  ■pltgr(tmax)' 
integer*4  tmax  ♦/ 

{ 


int  ij.tkval; 

setvideomode  (HRESBW); 
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_getvideoconfig  (5ivc); 
setcolor  (2); 
jnoveto  (52,172); 
lineto  (640,172); 
moveto  (52,172); 

_lineto  (52,0); 
jnoveto  (50,0); 

_lineto  (54,0); 
jnoveto  (50,43); 

^lineto  (54,43); 
moveto  (50,86); 

_lineto  (54,86); 

_(novetO  (50,129); 

Jineto  (54,129); 
tkval  =  (♦tmax+10)/10; 
for  (i=52;  i<=632;  i=i+58)  { 
jnoveto  ( 1,1 70); 

_lineto  (i,174); 

} 

_settextposition(l,3); 
printfCIOOO*); 
^settextposition(6,4); 
printfClOO’); 
jettextposition(11,5); 
printf ("10'); 
settextposition(17,6); 
printfCT); 
_settextposit1on(22,5); 
printfC.I*); 
_settextposition(23,7); 
pnntf(‘0“); 
settextposit1on(10,1); 
pnntfCHCr); 
jettextposition(  11,1); 
printfCppm*); 


for  (i  =  1;  i<=n);  i++  )  ( 

settextposition  (  23,(7*(i+1))  ); 
printf  (a‘,{tkva1+(i'1)*tkval)  ); 

} 

settextposition(25,35): 
printf  {‘time  (min)‘); 

) 

Sinclude  <stciio.h> 

Kinclude  <math.h> 

Sinc'lude  <qraph.h> 

struct  videoconfig  vc;  /*  variable  vc  of  type  videoconfig  */ 

/♦  This  function  accepts  the  elapsed  time  and  the  concentration  in 
ppm  and  plots  it  on  a  log  plot,  tmax  is  the  time  scale  (min) 
maximum  */ 

/*  This  C  funtion  plots  properly  in  the  _HRES8W  mode  and  should  be 
used  with  "pltgr“  */ 

/»  called  as  FORTRAN  subroutine  ■pltpts(etim,cppm,tffiax)* 
real*4  etim.cppm 
integer*4  tmax  */ 

void  pltpts(etim,cppm,tmax) 
int  near  *tmax; 
float  near  ♦etim; 
float  near  »cppm; 

{ 

int  xscale; 
int  x,y; 

xscale  =  10*((*tmax+10)/10); 

X  =  (int)  ((548. OT  *etim)/xscale); 
if(*cppm  <=  0.1) 
y  =  0; 

else  if(*cppm  >=  100.0) 
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y  =  172: 


else 

I  y  =  (int)  (43.0*1og10{*cppiii)+43.0); 

_setpixel(x+52, 172-y); 

) 

Jinclude  <stdio.h> 

I  Jinclude  <graph.h> 

I 

void  rscrnO 

/*  This  C  function  clears  the  screen  and  sets  the  video  mode  to 
I  TEXTBWSO  (the  normal  text  mode  for  the  Compaq)  */ 

/*  called  as  FORTRAN  subroutine  'rscrn*  ♦/ 

I 


_clearscreen( ^GCLEARSCREEN) ; 

_setvideomode(  JEXTBW80); 

) 

;  assembly  function  to  clear  the  screen,  should  word  in  any  video  mode 
;  called  as  FORTRAN  subroutine  'CIS' 


.model 

large 

.code 

PUBLIC 

CLS 

CIS 

PROC 

EXTRN 

HOME: FAR 

PUSH 

BP 

mov 

ax.OfOOh 

int 

OlOh 

mov 

bh,07 

cmp 

al,03 

jle 

bhset 

cmp 

al.07 

je 

bh^set 
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xor 

bh,bh 

bh^set;  mov 

ax,0600h 

xor 

cx,cx 

mov 

dx,184Fh 

int 

OlOh 

call 

HOME 

POP 

BP 

RET 

4 

CLS  ENOP 

END 

assembly  function  to  position  the  cursor  at  the  top  lefthand  corner 
of  the  screen 

called  by  assembly  function  'CLS' 


.model  large 

.code 

public 

HOME 

HOME 

proc 

push 

bp 

mov 

ax,0200h 

xor 

bx,bx 

xor 

dx,dx 

int 

OlOh 

pop 

bp 

ret 

HOME 

endp 

end 

:  assembly  function  to  fetch  the  date  for  FORTRAN  subroutine  'getdat(dc;te)' 
;  integer*4  date 

.model  large 
.code 

PUBLIC  gdate 
gdate  PROC 


S6 


PUSH 

MOV 


BP 


BP,SP 

xor  AH, AH 

mov  AH,2AH 

int  21H 

les  8X, DWORD  PTR  [8P+6] 

mov  £S:(8X],DX 

mov  ES:[8X]+2,CX 

POP  BP 

RET  4 

gdate  ENDP 
END 

;  assembly  function  to  fetch  the  time  of  day  for  FORTRAN  subroutine 
;  ‘gettim(time)' 

;  integer*4  time 


•MODEL  LARGE 
•  CODE 

PUBLIC  gtime 
gtime  PROC 
PUSH 
MOV 

xor 

mov 

int 

1es 

mov 

mov 

POP 

RET 

gtime  ENDP 


BP 

BP.SP 

AH,  AH 

AH,2CH 

21H 

BX, DWORD  PTR  [BP+6] 

ES:[8X],DX 

ES:[BX]+2,CX 

BP 

4 
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END 

:  assembly  function  to  Input  a  low  byte  from  a  port  called  as  FORTRAN 
;  subroutine  “inplbCport, value)' 

;  integer’4  port, value 

.model  large 
.code 

PUBLIC  INPLB 
INPLB  PROC 

PUSH  BP  ;Saves  framepointer  on  stack 
MOV  BP,SP 

LES  BX, DWORD  PTR[BP+10] 

MOV  0X,ES;tBX] 

LES  BX, DWORD  PTR[8P+6] 

IN  AL,DX 

XOR  AH, AH 

MOV  £S:[BX],AX 

POP  BP  ;  Restore  the  framepointer 

RET  8 

INPLB  ENDP 
END 

;  assembly  function  to  output  a  low  byte  to  a  port  called  as  FORTRAN 
:  subroutine  *outlb(port, value)' 

;  integer*!  port, value 

.model  large 
.code 

PUBLIC  OUTLB 
OUTLB  PROC 

PUSH  BP  ;Saves  framepointer  on  stack 
MOV  BP,SP 
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LES  BX, DWORD  PTRiBP+10] 

MOV  OX,ES:[BX] 

LES  BX, DWORD  PTR[BP+6] 

MOV  AL,ES;[BX] 

OUT  DX.AL 

POP  BP  •,  Restore  the  framepointer 

RET  8 

0UTL8  ENOP 
END 

;  assembly  function  to  output  a  character  string  to  the  current  text 
;  cursor  position.  String  must  be  terminated  with  a  character 
;  called  as  a  FORTRAN  subroutine  ■pstrng(c)‘ 

;  character  c*(up  to  80) 

.model  large 
.code 

PUBLIC  pstrng 
pstrng  PROC 

PUSH  BP 

MOV  8P,SP 

les  BX, DWORD  PTR  [BP+6] 

mov  DX,BX 

push  OS 

mov  AX,ES 

mov  OS, AX 

mov  AH, 9, 

int  21H 

pop  OS 

POP  BP 

RET  4 

pstrng  ENOP 

69 


END 

hcl.obj:  hcl.foP 


fl  /c  hcl.for 

hciscrn.obj:  hclscrn.for 

fl  /c  hclscrn.for 

crtcis.obj:  crtcls.asm 

masm  crtcis; 

crt^home.obj;  crt  home. asm 

masm  crt^home; 

pstring.obj;  pstring.asm 

masm  pstring; 

gettime.obj:  gettime.for 

fl  /c  gettime.for 

gtime.obj:  gtime.asm 

masm  gtime; 

twait.obj:  twait.for 

fl  /c  twait.for 

rdparm.obj:  rdparm.for 

fl  /c  rdparm.for 

stpclr.obj:  stpclr.for 

fl  /c  stpclr.for 

setad.obj:  setad.for 

fl  /c  setad.for 

rdatod.obj:  rdatod.for 
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fl  /c  rdatod.for 


srSIfiO.obj:  sr9600.for 

fl  /c  sp9600.for 

atodfUe.lib:  stpclr.obj  setad.obj  rdatod.obj  sr9600.obj 

lib  atodfile  -+stpclr.obJ  -+setad.obJ  -+pdatod.obJ  -+sp9B00.obJ; 

sorter. obj:  sorter. for 

fl  /c  sorter. for 

analOl.obj:  analOl.for 

fl  /c  analOl.for 

wait. obj:  wait. for 

fl  /c  wait. for 

pltfr.obj;  pltgr.c 

qcl  /c  /AL  /Od  /FPi87  pltgr.c 

pltpts.obj;  pltpts.c 

qcl  /c  /AL  /Od  /FPi87  pltpts.c 

rscreen.obj:  rscreen.c 

qcl  /c  /AL  /Od  /FPi87  rscreen.c 

keyhit.obj:  keyhit.c 

qcl  /c  /AL  /Od  /FPi87  keyhit.c 

cfiles.lib:  pltgr.obj  pltpts.obj  rscreen.obj  keyhit.obj 

lib  cfiles  -tpltgr.obj  -+pltpts.obj  -trscreen.obj  -+keyhit.obj; 

inplb.obj:  inplb.asm 

masm  inplb: 


71 


outlb.obj: 


out  lb. asm 


masm  outlb; 

afiles.lib:  cptcls.obj  crthome.obj  pstring.obj  gtime.obj  inplb.obj  outlb.obj 

lib  afiles  -+crt_cls.obj  -+crt  home.obj  -tpstring.obj  -+gtime  -+inplb.obj  -+outlb.obj: 

ffiles.lib:  hclscrn.obj  twait.obj  gettime.obj  rdparm.obj  sorter. obj  analOl.obj  wait.obj 

lib  ffiles  -+hclscrn  -+twait  -+gettiine  -+rdparm  -+sorter  -+anal01  -+wait; 

hcl.exe:  hcl.obj 

link  /NOE  hcl , , .ffiles+atodf ile+afiles+cf 1les+c:qUickciib^libc7.lib: 
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APPENDIX  C 


OPTICAL  TRAIN 

The  optical  train  of  the  HCl  detector  is  shown  In  Figure  C.l.  Light 
from  the  HCl  lamp  is  brought  to  focus  at  the  center  of  the  rotating  matched 
filter  gas  cell  using  a  50  mm  focal  length  calcium  fluoride  lens  and  a 
steering  mirror.  The  HCl  lamp  is  imaged  by  using  a  beam  stop  which  limits 
the  aperture  to  the  center  portion  of  the  lamp.  This  stop  Is  placed  at  the 
focus  of  the  rear  window  of  the  lamp,  thus,  limiting  the  off-axis  rays  from 
the  sidewalls.  This  beam  stop  acts  as  somewhat  of  an  aperture  stop  for 
rays  in  front  of  the  rear  windov;,  but  does  not  limit  the  light  from  the 
center  of  the  lamp. 

Light  passing  through  the  gas  cell  is  then  directed  by  a  second 
steering  mirror  through  a  38mn  f.l.  lens  to  focus  the  image  of  the  lamp 
into  the  White  cell.  The  White  cell  consists  of  three  mirrors,  each  having 
a  radius  of  curvature  of  250  mm  and  arranged  in  the  conventional  White  cell 
geometry.  The  center  of  curvature  of  the  two  rear  mirrors  have  been 


Figure  C.l.  Schematic  of  Optical  Layout. 
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machined  to  be  separated  by  3  mm.  This  separation  provides  six  spots,  each 
separated  by  6  mm  on  the  top  portion  of  the  front  mirror  before  the  beam 
exits  off  the  edge.  These  spots,  together  with  the  seven  spots  on  the 
bottom  portion  of  the  front  mirror  provide  a  total  of  28  passes  through  the 
cell  before  exiting  the  other  side.  Exiting  light  is  passed  through  a 
second  38  mm  f.l.  lens  so  that  the  circle  of  least  confusion  of  the 
focussed  light  is  centered  on  the  detector  by  the  final  steering  mirror. 

SYSTEM  ALIGNMENT 

The  system  is  prealigned  in  the  initial  assembly.  A  minimum  of 
realignment  should  be  necessary,  and  the  primary  service  should  involve 
cleaning  the  windows  and  mirrors  of  the  White  cell.  Access  to  the  White 
cell  is  obtained  by  removing  the  aerosol  heater  and  dust  cover,  then 
removing  the  8-32  socket  head  cap  screws  from  around  the  cell  cover,  and 
sliding  the  cell  cover  off.  The  mirrors  can  be  cleaned  with  methyl  alcohol 
and  a  clean  cotton  swab.  However,  should  a  mirror  become  damaged, 
realignment  must  accompany  replacement.  Since  the  entire  optical  train  is 
sealed  from  the  atmospheric  sampling  area,  the  primary  portion  requiring 
realignment  is  the  White  cell  mirror  assembly,  which  must  be  taken  out  and 
resurfaced.  Each  mirror  has  been  plated  with  a  nickel  coating,  polished 
and  then  gold  overcoated.  When  these  mirrors  are  severely  contaminated, 
they  should  be  removed  for  cleaning  and  resurfacing. 

The  White  cell  field  mirrors  can  be  removed  as  a  unit  by  removing  four 
6-32  cap  screws  bolted  to  the  side  panels  and  the  two  4-40  flat  head 
cap-screws  from  the  bottom.  The  entire  assembly  can  thus  be  cleaned  and 
mirrors  resurfaced.  Each  mirror  can  also  be  removed  Individually;  however, 
much  more  work  would  be  involved  to  realign  the  mirrors,  as  shims  have  been 
placed  for  final  alignment  of  the  White  cell  due  to  machining  tolerances  on 
the  brassboard  system.  These  should  be  saved  for  reassembly.  The  front 
mirror  can  be  accessed  by  removing  the  four  8-32  socket  head  cap  screws 
next  to  this  mirror.  This  removes  the  entire  White  cell  optical  bench. 
Three  6-32  screws  also  hold  this  mirror  in  place.  Good  nominal  alignment 


can  be  achieved  by  simply  reboltlng  each  subsystem  back  into  their  original 
positions.  However,  if  the  system  has  been  disassembled,  realignment  can 
be  achieved  as  follows: 

In  order  to  re-allgn  the  optical  train,  the  circulation  pump  must 
first  be  removed.  This  gives  access  to  the  optical  axis  of  the  system. 
Since  the  lamp  is  transparent,  a  small  HeNe  alignment  laser  is  easily 
placed  in  place  of  the  circulation  pump  to  provide  a  beam  of  collimated 
light  through  the  rear  end  of  the  lamp.  There  are  a  minimum  of  adjustments 
within  the  optical  train.  Placement  of  the  lenses  and  mirrors  have  been 
machined  to  nominal  positions.  Because  of  machining  tolerances  In  the 
brassboard  instrument;  however,  some  realignment  may  be  necessary  when  the 
White  cell  :*s  removed  and  replaced.  The  optical  axis  of  the  system  is  held 
at  2.00  inches  (50.8  mm)  above  the  optical  bench.  The  laser  should, 
therefore,  provide  a  light  beam  at  this  height  that  is  parallel  with  the 
optical  bench.  This  laser  beam  should  pass  through  the  left  hand  slot  of 
the  front  White  cell  mirror  and  fall  onto  the  center  of  the  left  rear 
mirror.  The  i..curn  from  this  rear  mirror  should  fall  directly  under  the 
inside  edge  of  the  right  hand  slot.  The  further  passes  through  the  White 
cell  should  resemble  these  shown  in  Figure  C.2.  Adjustments  of  the  White 
cell  consists  of  putting  shims  under  each  of  the  rear  mirrors  until  the  two 
rows  of  return  spots  are  parallel  with  the  optical  bench.  Slight 
adjustments  of  the  rear  mirrors  can  be  accomplished  by  moving  these  mirrors 
within  their  mounting  holes  before  tightening.  This  "slop"  would  be  keyed 
out  of  production  Instruments. 

Once  the  laser  beam  exits  from  the  right  hand  side  of  the  White  cell, 
it  should  strike  the  center  of  the  collection  lens  and  be  focused  onto  the 
detector.  Slight  misalignments  due  to  non-exact  placement  of  the  White 
cell  mirrors  may  be  compensated  for  by  adjusting  the  steering  mirror. 
Gross  adjustments,  however,  are  best  made  by  alignment  of  the  White  cell. 


FRONT  MIRROR 
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APPENDIX  D 


ELECTRONIC  DESIGN  FOR  THE  HCl  MONITOR 

The  electronic  circuits  for  opcratlnc  the  HCl  monitor  are  physically 
located  in  three  places: 

•  within  the  optical  module 

•  within  the  power  supply  module 

•  inside  the  portable  personal  computer 

The  circuits  are  briefly  described  in  the  following. 

A.  OPTICAL  MODULE  ELECTRONICS 

Figure  D.l  shows  the  circuits  used  to  amplify  and  buffer  the  detector 
signal,  to  control  the  detector  refrigeration  temperature,  to  regulate  the 
chopper  wheel  motor  speed,  and  to  regulate  the  HCl  lamp  temperature.  The 
aerosol  heater  operates  directly  from  a  5-volt  regulated  power  supply  and 
does  not  need  control  circuitry.  It  draws  a  current  of  about  5  amperes. 

All  the  circuits  in  the  optical  module  were  mounted  on  a  printed 
circuit  board  which  was  soldered  to  the  detector  as  an  integral  unit.  The 
power  supply  and  signal  paths  are  provided  by  an  18-pin  military  connector. 

1.  Uetfjctor  Amplifier  and  Buffers 

The  lead  selenlde  detector  operates  as  a  photo-conductive  device 
(Rjj)  with  a  resistance  of  3M  ohm.  To  obtain  maximum  signal  from  the 
detector,  the  load  resistor  Rg  is  chosen  to  match  this  Impedance.  The 
signal  is  AC  coupled  to  a  buffer  Ula  and  then  amplified  by  100  by  Ulb, 
and  buffered  again  by  U2a.  The  output  of  U2a  is  used  to  drive  the 
motor  controller  and  the  output  driver  circuit. 
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The  output  driver  (U2c,  U2d  and  Q2)  amplifies  the  signal  by  2  and 
adds  to  it  a  positive  offset.  The  latter  is  necessary  to  keep  the 
output  transistor  Q2  from  cut-off.  The  offset  adjustment  is  done  by 
trimming  RIO  so  that  the  base  line  of  the  output  signal  is  above 
cut-off  (~0.2v  above  ground).  The  present  circuit  can  drive  1000  ft 
of  RG-58-U  cable  with  a  50  ohm  termination  with  little  signal 
degradation . 

2.  Detector  Refrigerator  Controller 

The  lead  selenide  detector  is  cooled  internally  by  a 
thermo-electric  device.  The  temperature  of  the  detector  is  provided 
by  an  internal  thermister,  the  calibration  curve  of  which  is  shown  in 
Figure  D.2.  To  maintain  a  constant  temperature,  the  current  to  the 
thermo-electric  cooler  Rf  is  regulated  by  a  negative  feedback 
amplifier  composing  of  L)2b  and  Qj .  The  feedback  maintains  the 
thermister  resistance  at  lOKQ  which  corresponds  to  a  temperature  of 
-26°C.  The  circuit  has  been  tested  satisfactorily  over  an 
environmental  temperature  range  of  20  to  70°C. 

3.  Motor  Speed  Controller 

The  detector  signal  from  buffer  U2a  is  filtered  by  R30  and  C13 
and  AC  coupled  into  the  frequency  to  voltage  converter  U3 .  The 
converted  voltage  is  compared  with  a  reference  voltage  set  by  R19,  and 
the  difference  signal  is  used  in  a  feedback  loop  to  maintain  the  motor 
speed.  Different  motor  speeds  may  be  obtained  by  adjustment  of  R19. 
The  normal  speed  is  set  at  -1900  rpm. 
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Figure  D.2.  Calibration  Curve  of  Thermistor  in  the  PbSb  Detector. 
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4.  The  HCl  Lamp  Temperature  Controller 

The  temperature  of  the  HCl  lamp  is  measured  by  a  Type  K 
thermocouple  located  at  about  the  middle  of  the  lamp.  The 
thermocouple  voltaj’e  is  amplified  by  U4  (which  also  provides  an 
ice -point  internal  reference).  This  voltage  is  compared  with  a  preset 
value  set  by  R29.  The  difference  is  used  to  drive  the  feedback 
amplifier  U5a  and  to  provide  the  necessary  lamp  current  for  stable 
temperature  operation.  The  temperature  to  voltage  transfer  function 
for  114  is  given  in  the  following  table. 


1(C) 

V(Volts) 

500 

5.107 

600 

6.101 

700 

7.206 

(source;  Analog 

800 

8.232 

Devices  Data  Book) 

900 

9.233 

A  typical  lamp  temperature  history  is  shown  in  figure  I). 3.  The 
heat  up  process  started  at  A.  (The  lamp  was  slightly  hotter  than  room 
temperature  because  it  had  not  cooled  down  completely  from  a  previous 
test.)  The  heat  up  time  (to  B)  was  about  3  1/2  minutes  with  a  slight 
temperature  overshoot  at  B.  At  C,  the  temperature  set  point  was 
readjusted  (by  adjustment  of  K29),  and  the  lamp  temperature  increased 
to  U.  The  temperature  remained  constant  for  the  rest  of  the  test 
period . 

B.  POWKK  SUPPLY  MODULE 

The  power  supply  module  houses  five  regulated  power  supplies  (see 
Figure  D.4).  The  lOO-volt  low  current  supplv  for  the  detector  was 
regulated  by  a  zener  diode.  The  other  four  supplies  were  commercial  units. 
The  25-volt  supply  was  a  28  volt  commercial  unit  adjusted  down  to  25  volt. 
The  AC  line  current  was  about  2  amperes. 
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C .  COMPUTER  MODULE 


The  computer  module  consists  of  a  Compaq  portable  PC-XT  compatible 
personal  computer.  The  analog  to  digital  conversion  is  performed  by  a 
Data-Translation  DT  2812  board.  The  signal  from  the  HCl  monitor  is 
buffered  by  the  circuit  shown  in  Figure  D.5.  Also  shown  on  the  figure  is 
the  circuit  which  derives  a  clock  signal  from  the  analog  signal  to  strobe 
in  the  data. 

Referring  to  Figure  U.5,  the  analog  signal  is  filtered  by  R102  and 
C102,  buffered  by  UlOlb,  and  is  then  fed  to  the  analog  input  of  the  A  to  D 
converter.  The  clock  signal  is  derived  from  the  analog  signal  using  a 
Schmidt  trigger  (UlOlc),  the  output  of  which  is  inverted  by  QlOl  and 
connected  to  the  monostable  mul tlvibrators  U102  and  U103.  The  time  delay 
of  the  trigger  signals  (which  can  be  measured  at  pin  5  of  U102b  and  U103b) 
from  the  trailing  edge  of  the  signal  (at  the  inputs  (pin  1)  of  U102a  and 
U103a)  are  set  by  R109  and  R112  .respectively .  The  timing  diagram  is  shown 
at  the  bottom  of  Figure  D.5.  The  trigger  signals  are  combined  by  a  NOR 
gate  (U104)  and  inverted  by  Q102  to  provide  the  strobe  signal  for  the  A  to 
D  conversion. 

The  +12-,  -12-,  and  +5-volt  power  supplies  for  the  above  circuit  were 
obtained  from  the  bus  connector  of  the  computer. 

D.  ADJUSTMENTS 

The  adjustments  within  the  optical  module  are; 

•  Lamp  temperature  adjustment  R29 . 

•  Motor  speed  adjustment  R19. 

•  Signal  offset  adjustment  Rll. 

The  locations  of  these  adjustment  trim  pots  are  shown  in  Figure  D.6.  The 
adjustments  within  the  computer  module  are: 

•  Time  delay  adjustment  for  the  signal  strobe  R109. 

•  Time  delay  adjustment  for  the  baseline  strobe  RH2. 

The  location  of  these  adjustment  trim  pots  are  shown  in  Figure  D.7. 


84 


BOTTOM 


TOP 


Figure  D.6.  Adjustments  on  PC  Board  Within  the  Optical  Module. 


Figure  D.7.  Adjustments  on  PC  Board  Within  the  Computer  Module 
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